COEFFICIENTS  FOR  CALCULATING 
RADIATION  IMPEDANCES  AND  FAR-FIELD 
PRESSURES  OF  FREE-FLOODED  RING 
TRANDUCERS 

Peter  H.  Rogers,  et  al 

Naval  Research  Laboratory 
Washington,  D.  C. 

14  June  1974 


DISTRIBUTED  BY: 


National  Technical  Information  Service 
U.  S.  DEPARTMENT  OF  COMMERCE 

5285  Port  Royal  Road,  Springfield  Va.  22151 


SECURITY  CLASSIFICATION  OF  This  PAGE  flWiwi  Data  Entarad) 


REPORT  DOCUMENTATION  PAGE 


.  REPORT  NUMBER  2.  GOVT  ACCESSION  NO. 

NRL  Report  7749 


4.  TITLE  'and  Subtitle) 

COEFFICIENTS  FOR  CALCULATING  RADIATION 
IMPEDANCES  AND  FAR-FIELD  PRESSURES  OF  FREE- 
FLOODED  RING  TRANSDUCERS 


7.  AUTHORf*; 

Peter  H.  Rogers  and  Joseph  F.  Zalesak 


S.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Naval  Research  Laboratory 
Washington,  D.C.  20375 


II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Department  of  the  Navy 
Office  of  Naval  Research 
Arlington,  Virginia  22217 


MONITORING  AGENCY  NAME  A  ADDRESS  (It  dlllarant  from  Controlling  Olllca) 


■ 

I 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


1.  RECIPIENT’S  CATALOG  NUMBER 

AC-  7?2  35V 


5.  TYPE  OF  REPORT  4  PERIOD  COVERED 

Interim 


(.  PERFORMING  ORG.  REPORT  NUMBER 


(.  CONTRACT  OR  GRANT  NUMBERf*; 

S01-29 

S01-30 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  4  WORK  UNIT  NUMBERS 

RR1 31-03-41-5226 
RR1 31-03-41-5227 


U.  REPORT  DATE 

June  14, 1974 


13.  NUMBER  OF  PAGES 
128 


IS.  SECURITY  CLASS,  fof  thla  raport) 

Unclassified 


IS*.  DECLASSIFICATION/ DOWN  GRADING 
SCHEDULE 


IS.  DISTRIBUTION  STATEMENT  (ol  thla  Raport) 

Approved  for  public  release;  distribution  unlimited. 


<7.  DISTRIBUTION  STATEMENT  (ol  tha  ab, tract  antarad  In  Block  30,  II  dlllarant  from  Raport) 


It.  KEY  WORDS  ( Conttnu a  on  ravataa  aid a  li  nacaaaary  and  Identity  by  block  numbar) 

Acoustic  pressure  Free-flooded  rings  Piezoceramic  cylinders 

Acoustic  radiatioi  Impedance  (acoustical)  Sound  transducers 

Computer  programs  Influence  coefficients  Underwater  transducers 

Far-field  pressure  patterns  Omnidirectional  transducers 


20.  ABSTRACT  (Contlnua  on  ravataa  alda  H  nacaaaary  and  identity  by  block  ntimbor) 

Tables  of  coefficients  for  determining  the  radiation  impedance  and  far-field  pressure  of  a 
free-flooded  vibrating  ring  in  an  arbitrary  medium  at  45  frequencies  for  each  of  36  ring  geo¬ 
metries  were  produced.  The  only  restriction  on  the  known  velocity  distribution  for  using  the 
tables  is  that  the  inside,  top  and  bottom,  and  outside  normal  surface  velocities  must  be  uniform. 
Simple  formulas  for  determining  the  normal  surface  velocity  distribution  for  some  special  cases 
and  formulas  for  determining  the  radiation  resistance  and  far-field  pressure  in  the  low-frequency 
limit  were  developed. 


DD  1473  EOITION  OF  I  NOV  «5  IS  OBSOLETE 

S/N  0  102-014-  6601 

Rt'produco'l  l'y 

NATIONAL  TECHNICAl 
INFORMATION  SERVICE 

(I  Ofpnrtniont  of  Ommrrce 
Sonngfield  VA  22151 


SECURITY  CLASSIFICATION  OF  THIS  RACE  fBfc.n  Data  tntarad) 


CONTENTS 


INTRODUCTION .  1 

DESCRIPTION  AND  USE  OF  TABLES .  1 

LOW-FREQUENCY  VELOCITY  DISTRIBUTIONS .  5 

LOW-FREQUENCY  RADIATION  CHARACTERISTICS .  7 

CONCLUDING  REMARKS .  7 

REFERENCES  .  8 

APPENDIX  A  —  Impedance  and  Pressure  Coefficients .  9 

APPENDIX  B  —  NRL  SHIP  Program .  82 


m  Preceding  page  blank 


COEFFICIENTS  FOR  CALCULATING  RADIATION  IMPEDANCES  AND  FAR-FIELD 
PRESSURES  OF  FREE-FLOODED  RING  TRANSDUCERS 


INTRODUCTION 

The  free-field  radiation  characteristics  of  a  vibrating  free-flooded  ring  (right  rectan¬ 
gular  toroid)  are  completely  determined  by  the  frequency  of  oscillation,  the  ring  geom¬ 
etry,  the  normal  velocity  of  the  surface  of  the  ring,  and  the  density  and  sound  speed  of 
the  acoustic  medium.  The  parameters  of  the  medium  are  easily  factored  out  of  the  ex¬ 
pressions  characterizing  the  radiation  from  free-flood  rings.  Hence,  radiation  from  ring 
transducers  can  be  characterized  by  the  relative  dimensions  and  acoustic  size  of  the  ring 
independent  of  the  parameters  of  the  medium.  The  normal  velocity  distribution  on  the 
surface  of  a  ring  transducer,  however,  is  a  function  not  only  of  the  ring  geometry  and 
frequency  of  oscillation  but  also  of  the  ring  material  parameters,  the  type  of  drive,  and 
the  acoustic  load.  It  is  obviously  not  possible  to  include  all  these  parameters  in  a  table 
of  radiation  characteristics  for  free-flooded  rings.  Therefore,  as  a  compromise,  we  will 
present  numerical  data  corresponding  to  each  surface  having  a  unit  outward  velocity, 
with  the  remaining  surfaces  fixed  in  space.  Formulas  will  be  presented  for  combining  the 
data  to  obtain  numerical  values  corresponding  to  actual  physical  situations. 

The  data  presented  in  the  tables  in  Appendix  A  were  generated  with  the  latest  ver¬ 
sion  of  the  NRL  SHIP  program  [1],  a  fast  computer  algorithm  for  determining  the 
radiation  characteristics  of  free-flooded  ring  transducers  whose  normal  axisymmetric 
velocity  distribution  is  known.  A  listing  of  the  program  is  given  in  Appendix  B. 

Two  tables  are  presented  for  each  of  the  36  ring  geometries  considered.  The  first 
table  for  each  geometry  contains  impedance  coefficients,  and  the  second  table  contains 
pressure  coefficients.  The  coefficients  in  these  tables  are  unitless  complex  numbers.  The 
surface  velocity  distribution,  the  frequency  of  oscillation,  an  absolute  dimension  of  the 
ring,  and  the  density  and  sound  speed  of  the  medium  must  be  known  before  these  coef¬ 
ficients  can  be  converted  into  absolute  numbers.  In  the  case  of  the  far-field  pressure,  the 
distance  at  which  the  pressure  is  to  be  determined  must  also  be  known.  The  numbers  in 
the  tables  are  believed  to  be  accurate  within  1%,  except  for  those  that  determine  the  real 
part  of  the  radiation  impedance  for  very  thin  geometries  at  low  frequencies.  However,  at 
very  low  frequencies  the  radiation  resistance  (and  the  far-field  pressure)  can  be  determined 
from  simple  formulas. 


DESCRIPTION  AND  USE  OF  TABLES 

Figure  1  shows  the  ring  geometry  nomenclature  used  in  developing  the  tables.  T 
and  H  are  defined  by  the  identities 


Note:  Manuscript  submitted  March  27,  1974. 
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Fig.  1  — Ring  geometry  nomenclature 


T  =  —  (!) 

a 

and 

flni,  (2) 

where  t  =  thickness,  a  =  mean  radius,  and  h  =  ring  height.  Each  table  presents  impedance 
and  pressure  coefficients  for  a  specific  value  of  T  and  H  as  a  function  of  ka  where  k,  the 
wave  number,  is  defined  by 


k  =  M  (3) 

f  being  the  frequency  of  oscillation  and  c  the  velocity  of  sound.  The  36  geometries  con¬ 
sidered  correspond  to  all  possible  combinations  of  T  —  0.05,  0.1,  0.2,  0.3,  and  0.5  and 
H  =  0.1,  0.2,  0.5,  1.0,  2.0,  3.0,  and  5.0  as  well  as  the  special  case  T  =  2/9  and  H  =  4/9. 
For  each  geometry,  ka  =  0.01,  0.05,  0.10,  0.15,  0.20,  0.30,  0.40,  0.50,  ...,  3.70,  3.80, 
3.90,  4.00,  4.50,  and  5.00. 

The  subscripts  1,  2,  and  3  refer  to  the  inside,  top  and  bottom,  and  outside  surfaces, 
respectively.  The  positive  normal  direction  is  defined  to  be  outward  from  the  ring  mate¬ 
rial.  An  impedance  Z,-y  can  be  defined  as 


Zij 


PiAi 
Vi  pc  A 


(4) 


where  Pj  is  the  average  rms  pressure  on  surface  j  caused  by  a  normal  rms  velocity  Vi  on 
surface  i.  The  area  of  surface  j  is  Aj,  the  entire  surface  area  of  the  ring  is  A,  and  p  is  the 
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density  of  the  medium.  In  Eq.  (4)  it  is  assumed  that  the  velocity  Vj  is  constant  or  can 
be  approximated  by  a  constant  over  the  surface(s)  i.  The  total  complex  power  P  radiated 
is 


3 

P=  £  pcAViZijVf  (5) 

i 

where  the  asterisk  denotes  the  complex  conjugate. 

If  one  wishes  to  define  a  single  impedance  Z  for  the  entire  ring,  it  is  necessary  to 
choose  an  appropriate  reference  velocity  V.  For  example,  V  may  be  the  velocity  of  the 
mean  radius,  or  V\t  V2,  or  V3,  or  any  linear  combination  of  these.  Then,  by  definition, 

P  =  ZVV*  ,  (6) 

so  that 

3 

£  PcAViZijVf 


If 

Zl  s 
z2  = 
z3  * 

and 

then  the  radiation  impedance  is  given  by 

3 

£  Zi\Vi\2  +  Re  ( Z\  V2  V3*  +  Z'2 V3  V*  +  Z’3  Vx  V*)  . 
i-l  '  '  _ 

The  quantities  Z,  and  Z\  are  listed  in  the  tables  of  impedance  coefficients.  If  the  veloc¬ 
ities  are  real  and  the  reference  velocity  is  unity,  then  the  expression  for  the  radiation 
impedance  becomes 


pcA 

VV* 


vv* 

Zu  . 

Z22 . 

Z33 .  (8) 

^23  +  2 32  , 
z3l  +  Zi3  , 

Zl2  +  z2l  , 


Z  =  PcA(j2  Z‘Vi2  +  Z'lV2V*  +  Z2V3V1 


Z3V1 
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where  now  the  velocities  V,  are  unitless.  As  an  example,  suppose  we  have  an  isotropic 
ring  whose  dimensions  are  a  =  0.1  m,  h  =  0.2  m,  and  t  =  0.02  m,  and  we  want  to  know 
the  radiation  impedance  Z  in  water  at  2387  Hz.  For  this  example,  assume  that  the  rela¬ 
tive  surface  velocity  distribution  of  the  ring  is  that  of  the  ring  vibrating  in  a  vacuum  at 
its  first  radial  resonance.  Assume  further  that  this  velocity  distribution  can  be  approxi¬ 
mated  with  sufficient  accuracy  by 


Vi  =  -1  -  \  To, 

V2  =  -jHo, 
and 

V3  =  1  -  \  To 

where  a  is  Poisson’s  ratio  and  the  velocity  distribution  represents  surface  velocities  rela¬ 
tive  to  the  reference  mean  radial  velocity  [2].  If  a  is  chosen  to  be  0.3,  the  relative 
velocities  (unitless  numbers)  become  V\  =  -1.030,  V2  =  -0.300,  and  V3  =  0.970.  The 
impedance  coefficient  tables  in  Appendix  A  are  listed  according  to  the  values  of  T  and 

H.  For  each  geometry,  the  corresponding  table  lists  impedance  coefficients  as  a  function 
of  ka.  In  this  example,  T  =  0.2,  H  =  2.0,  and  ka  =  1.0.  The  problem  is  to  evaluate 

z  -  ?ca(£  z,v,2  +  z;v2v3  ♦  zjv.v,  ♦  zjv.v.j. 

For  water,  pc  is  1.5  X  106  kg/s  m2.  The  total  surface  area  A  of  the  ring  is  0.2765  m2. 
The  values  of  V)  have  been  determined  previously  and  the  impedance  coefficients  are 
listed  in  Table  A23a.  The  radiation  impedance  in  this  example  is  (1.627  X  106  +  /' 

I. 064 X  106)  kg/s  or  mks  mechanical  ohms.  The  convention  here  is  that  massive  react¬ 
ances  are  positive.  (Large  negative  reactances  can  occur  in  the  region  of  a  strong  cavity 
resonance;  i.e.,  the  radiation  reactance  becomes  springlike.) 

Interpolation  between  frequencies  listed  is  valid  as  long  as  the  impedances  are  not 
changing  too  rapidly.  Interpolation  between  geometries  can  be  used  to  indicate  trends. 
The  numerical  value  of  the  impedance  obtained  in  this  way  is  unlikely  to  be  accurate. 

The  pressure  coefficients  p,°  and  pf°  are  uritless  complex  numbers  from  which  the 
far-field  pressure  along  the  radial  and  axial  directions,  respectively,  can  be  calculated  ac¬ 
cording  to 


i  =  l 


where  d  is  the  distance  at  which  the  pressure  is  to  be  determined  and  the  values  of  V, 
are  the  actual  (not  relative)  surface  velocities  of  the  ring.  For  p°,  substitute  the  p,° 
coefficients;  for  p90,  substitute  the  pf°  coefficients.  The  only  pressure  coefficients 
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listed  in  the  tables  in  Appendix  A  are  those  corresponding  to  0°  and  90°;  however,  pres¬ 
sure  coefficients  have  been  determined  at  5  intervals  and  are  stored  on  magnetic  tape. 
(These  coefficients  are  available  on  request.)  The  reason  the  formula  for  far-field  pressure 
is  simple  is  that  linear  superposition  holds. 

In  the  previous  example,  assume  that  the  far-field  pressure  at  1  m  is  to  be  deter¬ 
mined  when  the  mean  radial  velocity  is  1  m/s.  In  this  case  the  relative  surface  velocities 
listed  are  the  actual  surface  velocities  and  the  quantity  pca/d  is  1.5  X  105.  The  far-field 
pressure  at  1  m  is  |p°|  =  6.192  X  105  Pa,  or  235.8  dB  re  1  pPa  at  1  m,  and  |p90|  = 

7.419  X  104  Pa,  or  217.4  dB  re  1  pPa  at  1  m. 

These  pressures  are  for  a  ring  vibrating  with  a  mean  radial  velocity  of  1  m/s,  an  ab¬ 
normally  large  value  for  the  usual  transducer  materials  and  a  ring  of  the  geometry  of  this 
example.  This  is  the  reason  for  the  high  far-field  pressure.  The  actual  pressure  decreases 
linearly  as  the  mean  radial  velocity  of  the  ring  is  reduced. 

The  values  of  the  resistive  portion  of  the  impedance  coefficients  at  low  frequencies 
for  rings  of  small  t/a  ratio  have  errors  somewhat  larger  than  1%.  However,  these  occur 
well  into  the  region  where  the  radiation  is  monopole.  A  simple  formula  for  predicting 
the  radiation  resistance  will  be  presented  in  the  next  section.  Also,  for  certain  combina¬ 
tions  of  ring  geometry  and  velocity  distribution,  a  large  subtraction  error  results  from 
using  the  tables  for  calculating  the  radiation  resistance  in  the  monopole  region.  Again 
the  cure  is  to  use  the  monopole  radiation  formulas  presented  in  the  next  section.  The 
tables  can  easily  be  extended  to  even  lower  frequencies  because  in  the  monopole  region 
the  radiation  reactance  and  the  magnitude  of  the  far-field  pressure  are  directly  propor¬ 
tional  to  the  frequency  and  the  radiation  resistance  is  proportional  to  the  square  of  the 
frequency.  The  same  guidelines  apply  to  pressure  interpolation  as  apply  to  impedance 
interpolation. 


LOW-FREQUENCY  VELOCITY  DISTRIBUTIONS 

At  frequencies  somewhat  below  the  first  mechanical  resonance  of  the  free-flooded 
ring,  the  ring  is  essentially  in  static  equilibrium  with  whatever  driving  force  is  causing  the 
ring  to  vibrate.  As  a  specific  example  we  will  consider  a  piezoceramic  ring  driven  in  three 
different  ways  [3,4].  The  first  case  is  a  ring  driven  through  its  thickness  t;  i.e.,  a  ring 
polarized  and  driven  with  electrodes  on  the  inside  and  outside  surfaces.  Only  one  con¬ 
stitutive  relation  is  needed  to  establish  the  velocity  distribution: 

S  =  s'  CJ  +  d(6 

where  S  is  the  strain,  sf"  is  the  compliance  matrix  at  constant  electric  field,  3  is  the 
stress,  d  is  a  matrix  giving  the  strain  per  applied  electric  field,  and  6  is  the  applied  elec¬ 
tric  field  [5].  For  low-frequency  operation,  assume  stat*c  equilibrium  at  any  instant  of 
time;  i.e.,  set  J  =  0.  Note  that  this  condition  does  not  preclude  the  possibility  of  a  pre¬ 
stressed  ring.  Thus, 


S  =  df£ . 
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The  usual  convention  is  to  let  the  3  direction  be  the  polarization  direction,  and  in  our 
case  it  will  also  correspond  to  the  direction  of  the  applied  field  6.  Since  at  any  fixed 
frequency  the  velocity  is  proportional  to  the  strain,  it  is  easily  shown  that  the  relative 
velocity  distribution  is  given  by 


and 


where 


V!  =  -1  -  \  T  -i,  , 
V2  -  \  »• 

Va  =  l  -  i  T  A 

£  a 

,  ^31 

a  =  — • 


The  second  case  to  consider  is  the  one  in  which  the  ring  is  driven  through  its  height 
h;  i.e.,  the  ring  is  polarized  and  driven  with  electrodes  on  the  top  and  bottom  surfaces. 
The  analysis  is  again  straightforward,  and  the  resulting  relative  velocity  distribution  is 


and 


where  again 


Vi  =  -1  +yT, 

V2  =  -iff-,* 

*  a 
V3  =  1  +  \  T 

,  <*31 

°  d33’ 


In  the  third  case  the  ring  is  circumferentially  driven.  This  can  be  achieved  if  the 
ring  is  divided  into  an  even  number  of  segments  separated  by  electrodes.  The  assembled 
ring  is  then  polarized  and  driven  with  these  electrodes,  resulting  in  circumferential  fields 
and  hence  circumferential  drive.  The  analysis  is  again  straightforward,  and  the  resulting 
velocity  distribution  is 


Vi  =  -1  -  y  To' , 
V2  =  -jHo', 
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V3  =  1  -jTo' 

where  again 

,  <*3l 

°  =  <*33  ' 

Another  important  case  to  consider  is  the  velocity  distribution  at  the  first  radial 
resonance  of  an  isotropic  ring.  If  it  were  possible  to  excite  such  a  ring  and  the  ring  were 
not  loaded  heavily  acoustically,  then  its  velocity  distribution  would  be  approximated  [2] 
by 


Vi  =  -1  -  \Tuy 
V2  =  -~Ha, 

and 

V3  =  1  -  yTo 

where  a  is  Poisson’s  ratio  for  the  material.  These  velocity  distributions  are  all  relative  to 
a  unit  mean  radial  velocity. 


LOW-FREQUENCY  RADIATION  CHARACTERISTICS 

At  sufficiently  low  frequencies,  the  radiation  from  free-flooded  rings  is  monopole. 
In  this  case  simple  formulas  can  be  developed  to  predict  the  radiation  resistance  and  the 
far-field  pressure.  The  radiation  resistance  R  is  given  by 

1  l2 

2 V2T  +  (Vi  +  V3)H  +  \  (V3  -  V\)HT 


R  =  pcira2(ka)2 


where  V\,  V2,  and  V3  are  relative  velocities;  i.e.,  unitless  relative  to  a  unit  mean  radial 
velocity.  For  a  mean  radial  velocity  V,  the  far-field  pressure  at  distance  d  is  given  by 


CONCLUDING  REMARKS 

The  formulas  presented  were  experimentally  tested,  and  agreement  was  quite  good 
[3,4].  A  novel  application  of  the  data  in  the  tables  is  to  measure  the  low-frequency 
directivity  patterns  as  a  function  of  frequency  for  a  piezoceramic  ring.  The  best  value  of 
o'  =  ~d3i/d33  that  accurately  predicts  the  quantity  |p°/p90|  is  a  good  measure  of  the 
ratio  of  the  two  parameters  d3 1  and  d33.  This  was  tried  and  the  agreement  with  pub¬ 
lished  values  was  again  quite  good  [3,4]. 
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Appendix  A 

IMPEDANCE  AND  PRESSURE  COEFFICIENTS 


LIST  OF  TABLES 

Table  Ala.  Impedance  Coefficients  for  T  =  0.05;//  =  0.1 
Table  Alb.  Pressure  Coefficients  for  T  =  0.05;  //  =  0.1 

Table  A2a.  Impedance  Coefficients  for  T  -  0.1;  //  =  0.1 
Table  A2b.  Pressure  Coefficients  for  T-  0.1;//  =  0.1 

Table  A3a.  Impedance  Coefficients  for  T  =  0.2;  H  =  0.1 
Table  A3b.  Pressure  Coefficients  for  T  =  0.2;  H  =  0. 1 

Table  A4a.  Impedance  Coefficients  for  T  =  0.3;  II  =  0.1 
Table  A4b.  Pressure  Coefficients  for  T  =  0.3;  H  =  0.1 

Table  A5a.  Impedance  Coefficients  for  T  =  0.5;  H  =  0.1 
Table  A5b.  Pressure  Coefficients  for  T  =  0.5;  H  =  0. 1 

Table  A6a.  Impedance  Coefficients  for  T  =  0.05;  H  ~  0.2 
Table  A6b.  Pressure  Coefficients  for  T  =  0.05;  H  =  0.2 

Table  A7a.  Impedance  Coefficients  for  T  =  0.1;//  =  0.2 
Table  A7b.  Pressure  Coefficients  for  T  =  0. 1 ;  H  =  0. 2 

Table  A8a.  Impedance  Coefficients  for  T  =  0.2;  H  =  0.2 
Table  A8b.  Pressure  Coefficients  for  T  =  0. 2;  H  =  0. 2 

Table  A9a.  Impedance  Coefficients  for  T  =  0.3;  H  =  0.2 
Table  A9b.  Pressure  Coefficients  for  T  =  0.3;  H  =  0.2 

Table  AlOa.  Impedance  Coefficients  for  T  =  0.5;//  =  0.2 
Table  AlOb.  Pressure  Coefficients  for  T  =  0.5;  H  =  0.2 

Table  A1  la.  Impedance  Coefficients  for  T  =  0.05;  H  =  0.5 
Table  Allb.  Pressure  Coefficients  for  T  =  0.05;//  =  0.5 

Table  A12a.  Impedance  Coefficients  for  T  =  0.1;  H  =  0.5 
Table  A12b.  Pressure  Coefficients  for  T  =  0.1;  H  =  0.5 

Table  A  13a.  Impedance  Coefficients  for  T  =  0.2;  H  =  0.5 
Table  A13b.  Pressure  Coefficients  for  T  =  0.2;  H  =  0.5 

Table  A14a.  Impedance  Coefficients  for  T  =  0.3;  H  =  0.5 
Table  A14b.  Pressure  Coefficients  for  T  =  0.3; //  =  0.5 

Table  A15a.  Impedance  Coefficients  for  T  =  0.5;  //  =  0.5 
Table  A15b.  Pressure  Coefficients  for  T=  0.5;//  =  0.5 

Table  A16a.  Impedance  Coefficients  for  T=  0.05;//  =  1.0 
Table  A16b.  Pressure  Coefficients  for  T  =  0.05;  H  =  1.0 

Table  A17a.  Impedance  Coefficients  for  T  =  0.1;//  =  1  0 
Table  A  17b.  Pressure  Coefficients  for  T-  0.1;  H  =  1.0 

Table  A18a.  Impedance  Coefficients  for  T=  0.2;  H  =  1.0 
Table  A18b.  Pressure  Coefficients  for  T  =  0.2;  H  =  1.0 


Table  A19a.  Impedance  Coefficients  for  T  =  0.3;  //  =  1.0 
Table  A19b.  Pressure  Coefficients  for  T  =  0.3;  //  1.0 

Table  A20a.  Impedance  Coefficients  for  T  =  0.5;  H  =  1 .0 
Table  A20b.  Pressure  Coefficients  for  T  =  0.5;  H  =  1.0 

Table  A21a.  Impedance  Coefficients  for  T  =  0.05;  H  -  2.0 
Table  A21b.  Pressure  Coefficients  for  T  =  0.05;  //  =  2.0 

Table  A22a.  Impedance  Coefficients  for  T  =  0.1;  H  =  2.0 
Table  A22b.  Pressure  Coefficients  for  T  =  0.1;  H  =  2.0 

Table  A23a.  Impedance  Coefficients  for  T  =  0.2;  H  =  2.0 
Table  A23b.  Pressure  Coefficients  for  7’  =  0.2;  H  =  2.0 

Table  A24a.  Impedance  Coefficients  for  T  =  0.3;  H  =  2.0 
Table  A24b.  Pressure  Coefficients  for  T  =  0.3;  H  =  2.0 

Table  A25a.  Impedance  Coefficients  for  T  =  0.5;  H  ~  2.0 
Table  A!;5b.  Pressure  Coefficients  for  T  =  0.5;  H  =  2.0 

Table  A26a.  Impedance  Coefficients  for  T  =  0.05;  //  =  3.0 
Table  A26b.  Pressure  Coefficients  for  T  -  0.05;  //  =  3.0 

Table  A27a.  Impedance  Coefficients  for  T-  0.1;//=  3.0 
Table  A27b.  Pressure  Coefficients  for  T  =  0.1;//  -  3.0 

Table  A28a.  Impedance  Coefficients  for  T  =  0.2;  H  =  3.0 
Table  A28b.  Pressure  Coefficients  for  T  =  0.2;  H  =  3.0 

Table  A29a.  Impedance  Coefficients  for  7’  =  0.3;  H  =  3.0 
Table  A29b.  Pressure  Coefficients  for  T  =  0.3;  H  =  3.0 

Table  A30a.  Impedance  Coefficients  for  T  =  0.5;  H  =  3.0 
Table  A30b.  Pressure  Coefficients  for  T  =  0.5;  H  =  3.0 

Table  A31a.  Impedance  Coefficients  for  T  =  0.05;  II  =  5.0 
Table  A31b.  Pressure  Coefficients  for  T  =  0.05;  H  =  5.0 

Table  A32a.  Impedance  Coefficients  for  T  =  0.1 ;  H  =  5.0 
Table  A32b.  Pressure  Coefficients  for  T  -  0.1;  H  =  5.0 

Table  A33a.  Impedance  Coefficients  for  T=  0.2;  //  =  5.0 
Table  A33b.  Pressure  Coefficients  for  T  =  0.2;  H  =  5.0 

Table  A34a.  Impedance  Coefficients  for  T  0.3;  H  =  5.0 
Table  A34b.  Pressure  Coefficients  for  T  =  0.3;  //  =  5.0 

Table  A35a.  Impedance  Coefficients  for  T  =  0.5;  H  -  5.0 
Table  A 35b.  Pressure  Coefficients  for  T-  0.5;  //  5.0 

Table  A36a.  Impedance  Coefficients  for  T-  2/9;  //  4/9 
Table  A36b.  Pressure  Coefficients  for  T  =  2/9;  H  =  4/9 
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Appendix  B 

NRL  SHIP  PROGRAM 


The  version  of  the  SHIP  program  used  to  generate  the  data  in  the  tables  differs  from 
the  version  of  the  SHIP  program  listed  in  NRL  Report  7240*  in  three  significant  ways. 
The  first  difference  is  the  inclusion  of  Subroutines  CFINT,  EXI,  and  EXI2.  These  sub¬ 
routines  calculate  asymptotic  expressions  for  the  parts  of  the  integrals  truncated  in  Sub¬ 
routines  CSTM,  CTSM,  and  CSSM,  thus  giving  improved  accuracy.  However,  this  is  a 
time-consuming  process.  To  offset  this  increase  in  computer  execution  time,  most  of  the 
trigonometric  functions  needed  repeatedly  in  the  program  are  calculated  and  prestored  in 
Subroutine  CALTRIG.  With  this  addition,  the  program  runs  in  about  the  same  time  as 
before  Subroutines  CFINT,  EXI,  and  EXI2  were  added,  but  more  memory  space  is  re¬ 
quired.  The  third  major  change  in  the  program  is  the  addition  of  Subroutine  SOLZIJ. 
This  subroutine  calculates  impedance  coefficients  Z,  and  Z\  for  ring  transducers  (see 
Eq.  (8)). 

A  convenient  feature  added  to  the  program  is  the  inclusion  of  entry  point  FIELD  in 
Subroutine  SOLRING.  For  a  given  frequency  and  ring  geometry  it  is  not  necessary  to 
execute  the  entire  program  for  each  velocity  distribution.  A  single  call  to  SOLRING  (or 
SOLZIJ)  prestores  on  input/output  unit  10  everything  that  is  needed  for  any  velocity 
distribution  for  a  given  frequency  and  ring  geometry.  For  alternate  velocity  distributions, 
call  FIELD  to  avoid  recalculating  these  quantities. 

Numerous  minor  changes  were  made  in  almost  all  subroutines;  hence,  the  entire  pro¬ 
gram  is  listed  in  this  appendix.  Included  also  is  a  sample  output  generated  by  the  listed 
program  corresponding  to  the  example  listed  in  the  text.  This  illustrates  the  use  of 
SOLZIJ  and  FIELD. 


*P.  H.  Rogers,  “SHIP  (Simplified-Helmholtz-Integral  Program),  A  Fast  Computer  Program  for  Calculating 
the  Acoustic  Radiation  and  Radiation  Impedance  lor  Free-Flooded-Ring  and  Finite-Circular-Cylinder 
Sources,”  NRL  Report  7240,  June  19,  i 97 2. 
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NRL  REPORT  7749 


PROGRAM  SHIP 

BANK (0) ./6/./S  IORSCST/ 

COMMON/TOY/NOOI < ISYMt I  COR/P  I T/JmAX* JMAXH* ImAX/HLK 1/H*A*EK*P1 
COMMON/VELO/VEL/DIST/OIST/RAO/RIN.ROUT 
COMMON/MXD/JMXT.MXD/RRCC/RHOC.NPTS 
DIMENSION  VEL (60) « JMXT (60) 

TYPE  COMPLEX  VEL 

C  ROUT  IS  THE  OUTER  RADIUS  OF  THE  RING  R IN  THE  INNER  RADIUS 

C  EK  IS  THE  WAVE  NUMBER 

C  IMAX  IS  THE  NUMBER  OE  BANDS  ON  THE  TOP  AND  BOTTOM  OE  THE  RING 

C  JMAXH  IS  HALE  ThE  NUMBER  OE  HANDS  ON  THE  INSIDE  AND  OUTSIDE  SUREACES 

C  JMAXh  AND  IMAX  MUST  BOTH  BE  LESS  THAN  UR  EUUAL  TO  TEN 

C  NQO 1  SPEC  1 E IES  GAUSSIAN  QUADRATURE  ORDER  -  NUD1  MAY  BE  10.  20  OR  32 

C  NOD  1  =  32  UNLESS  CHANGED  IN  SHIP 

C  E OR  CYLINDER  SE i  RIN  =  0.  AND  ROUT  EQUAL  TO  THE  RADIUS 

C  H  IS  THE  HALE  HEIGHT  OE  THE  RING 

C  RHOC  IS  THE  PROOUCT  OE  THE  DENSITY  AND  THE  SOUND  VELOCITY  OF  THE  MEDIUM 

C  VEL ( I )  SPECIFIES  THE  (COMPLEX)  VELOCITY  ON  THE  1TH  BAND 

C  I S YM  =  ♦  1  =  SYMMETRY  ABOUT  Z  =  0«  -1  =  ANTISYMMETRY  AND  ISYM=0  =  NO  SYMMETRY 

C  FOR  Z(I.J)  CALL  SOLZIJ 

C  FOR  A  FIXED  FREQUENCY  AND  GEOMETRY  CALL  SOLRING  FOR  THE  FIRST 
C  VELOCITY  DISTRIBUTION  THEN  CALL  FIELD  FOR  ANY  SUBSEQUENT  VELOCITY 
C  DISTRIBUTIONS.  THIS  SAVES  CONSIDERABLY  IN  TIME. 

C  IE  BOTH  IMPEDANCE  COEFFICIENTS  AND  RADIATION  IMPEDANCE  ARE  DESIRED  CALL 

C  SOLZIJ  FIRST  THEN  CALL  FIELD  REPEATEDLY 

CALL  TIME 
RHOC  =  1.5E6 
I SYM  =  1 

JMAXH  =  10  *  IMAX  =  10 
FK  =  10.0 
RIN  =  0.09 
ROUT  =  0.11 
H  =  0.1 
CALL  SOLZIJ 
CALL  TIME 
VIN  =  -1.03 
VTOP  =  -0.3 
VOUT  =  0.97 
DO  100  I  =  1,  10 
VEL ( I )  =  VIN 
VEL (1*101  =  VTOP 
100  VEL (1*20)  =  VOUT 
CALL  FIELD 
CALL  TIME 

NPTS  IS  THE  NUMBER  OF  EAREIELD  POINTS  CALCULATED  BETWEEN  0  AND  90  DEGREES 
DIST  IS  THE  FARFIELD  DISTANCE  (1C0R=1)  IE  ICOR=0  FIELD  CALC.  AT  INFINITY 
ICOR  =  0 
NPTS  =  IB 
CALL  EARFLO 
CALL  TIME 
END 
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SUBROUTINE  S0LZ1J 


COMMON/TOY/NQD1.ISYM,1COR/RAO/RIN.ROU1/LC/LCMAX*LCMAXH/NOPH/NOPR 

COMMON/BLKl/H,A.FK,PI/PQR/POR/VELO/VEL/AREA/AREA 

COMMON/P  I  T/JHAX  •  JMAXH*  1MAX 

DIMENSION  VEL (60 ) *  AREA  <601*  ZIJI3.3).  PUW(30).  Z(3>«  VTEMP<60> 
TYPE  COMPLEX  VEL*  ZIJ*  POP*  Z*  VTEMP 
IF (WIN. EO. 0.0)  GO  TO  200 
I SYMTEMP  =  ISYM 
JM AX  =  JMAXH  ♦  JMAXH 
LCMAXH  =  JMAX  ♦  1M AX 
LCMAX  =  LCMAXH  ♦  LCMAXH 
DO  299  I  =  1*  LCMAX 
299  VTEMP ( I )  =  VEL(I) 

NOPR  =  1 
ISYM  -  1 

AREAT  =  (WIN  ♦  ROUT)  *  (H  ♦  H  ♦  ROUT  -  RIN) 

PRINT  TOO 

700  FORMAT ( 1H1 ) 

PRINT  701  .WIN. ROUT  .H.FK.NQD1. IMAX. JMAX  ?  3,  nqu1=*I3 

701  FORMAT  I4*  R IN  =  *F 7 . 3*  ROUT=»F7.3«  H=*F7.3*  TK-*F7.3  NUUl 


i#  IMAX=*I3«  JMAX=* 13//) 
DO  300  I  =  l*  3 
00  300  J  =  1*  3 
300  ZIJ(I*J>  =  0.0 
N 1  =  JMAXH  ♦  1 
N2  =  JMAXH  ♦  IMAX 
N3  =  N2  ♦  1 
DO  30S  1=1*  JMAXH 

305  VEL(I)  =  1*0 

DO  310  I  =  Nl.  N2 
310  VEL ( 1 )  =  0.0 

DO  306  1  =  N3*  LCMAXH 

306  VEL ( I )  =  0.0 

DU  311  N  =  1*  3 
IF (N.EQ.  1  >  CALL  SOLRING 
IF (N.EQ.2)  400*  401 

400  CONTINUE 

DO  350  1=1*  JMAXH 

350  VEL(I)  =  0.0 

DO  351  I  =  Nl.  N2 

351  VEL ( I )  =  1*0 
CALL  FIELD 

401  IF (N.EQ.3)  402*  403 

402  CONTINUE 

DO  352  I  =  Nl.  N2 


352  VEL ( l )  =  0.0 
DO  353  I  =  N3. 

353  VEL ( I )  =  1*0 
CALL  FIELD 

403  continue 
DO  315  I 
315  ZI J< 1 *Ni 
00  320  I 
320  Z  I  J  <  2  *N ) 

DO  325  I 


LCMAXH 


1,  JMAXH 
ZIJ(1*N)  ♦ 
Nl.  N  2 
/ I J  <2  *N)  ♦ 
N3.  LCMAXH 


POR  (  I  ) 
PUR  ( 1  ) 


*  AREA (  I ) 
«  AREA  (  I ) 
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325  ZIJOtN)  =  ZIJ(3.NI  ♦  PQR  ( I  >  *  ARt  A  <  I ) 

311  CONTINUE 

00  330  1  =  1*  3 
00  330  J  =  It  3 

330  ZlJ(ltJ)  =  ZU(ItJ)  •  2.0  *  CMPLX(0.0t  -FK>  /  AREAT 
Zll)  =  ZlJ(2t3>  ♦  ZIJ(3t2l 

Z (21  =  ZlJ<3tl>  ♦  ZlJ(lt3) 

Z (3)  =  ZIJ(lt2)  ♦  ZIJ(2.1> 

PRINT  99 

99  FORMAT ( lXt*TH£  IMPEDANCE  COEFFICIENTS  FOR  THIS  RING  TRANSDUCER  IN 
JUNITS  OF  RHO  C  A  ARE* t / 1 

PRINT  lOOt  ZlJUtllt  ZIJ(2t2)t  ZIJ(3t3)t  Z(l)t  Z  (2  >  *  Z(3) 

100  FORMAT (2X.SHZ1  =  tC (E 14 .5 tE 1 A. 5) */ »2X 1 5HZ2  =  tC (E l4.5tEl4.5) t/t 

$  2Xt5HZ3  =  tC(E14.5tE14.5) t/tlXt6HZlP  =  tC <E 14 .StE14.5) t/ t 

5  lXtfeHZ2P  =  tC(E14.5.E14.5) ,/,lX,6HZ3P  =  t C <E 14.5tE14.5> t / ) 

NOPR  =  0 

ISYM  =  ISYMTEMP 

00  298  I  =  It  LCMAX 
298  VEL(I)  =  VTEMP ( 1 ) 

GO  TO  210 
200  PRINT  700 

101  FORMAT (lXt*THIS  PROGRAM  OOES  NOT  COMPUTE  IMPEDANCE  COEFFICIENTS  FO 
$R  CYLINDER  TRANSDUCERS*) 

210  CONTINUE 
RETURN 
END 


85 


ROGERS  AND  ZALESAK 


SUBROUTINE  SOLRING 

TYPE  COMPLEX  V » SUM J.RHOC, ANS, VEL  , l)M,GM  »  VtL  1  ,  ANSI  ,  ANS2, PUR 
DIMENSION  REST (6000) . ANS <60 ) . V < 30 >  » VEL (60) ,GM (30,60) *DM<bO*feO> 

1 .VEL1 < 30) ,PZ<20>  ,PLZ (20) , PR (20) .PLR120)  .ANSI  < 30 ) . ANS2 < 30 )  .AREA (60) 
DIMENSION  JMX  T (60 ) .  PQR(30) 

COMMON/LC/LCMAX .LCMAXH/BLKA/PR.PLW.PZ , PLZ/RCC/ JTOP 
C0MM0N/T0Y/NQU1 . ISYM.NQD3/PI T/JMAX. JMAXH. IMAX 
COMMON /HLK 1/H.A.EK.P 1 /DEL/OELH .OELZ/RRCC/RRCC  » NPTS 
COMMON/HAD/R IN.ROUT/MXD/JMXT .MXU/S/DM/ ANS/ ANS/ VELO/VEL 
C0MM0N/6/GM.REST/T 1DY/EAST/EPS/EPS/EPS1/EPS1/J0E/JT0P1 
COMMON/AREA/AREA/nOPK/NOPR/PQR/PQR 

DATAfNODl  =  32), (Pi  =  3. 1415926S359)  ,  (EPS  =  0.001), (EPS1  =  0.0001) 
DATA ( JT0P1  =  10)  ,  ( JTOP  =  10) 

JMAX  =  2*JMAXH 

EAST  =  PI  •  NODI  /  (8.0  *  EK  •  (ROUT  ♦  PI  •  H>> 

IE (RIN.EQ.O.O)  GO  TO  200 
LCMAXh  =  JMAX  ♦  IMAX 
A  =  (WIN  ♦  ROUT)  •  O.S 
GO  TO  201 

200  LCMAXH  =  JMAXh  ♦  IMAX 
A  =  ROUT 

201  LCMAX  =  LCMAXH  ♦  LCMAXH 
CALL  GQC 

CALL  ERING 

RHOC  =  -EK  *  CMPLX (0.0,  RRCC ) 

JO  14  J=l, LCMAXH 
11  ANSI ( J ) =  ANS2 ( J)  =  <0.0, 0.0) 

DO  14  1  =  1. LCMAXH 

DM ( I , J )  =  OM(I.J)  »  DM  <  I  ,  LCMAX  ♦  1  -J) 

GM(I,J)  =  ('iM(I.J)  ♦  GM  <  I ,  LCMAX  ♦  1  -J) 

14  CONTINUE 
REWIND  10 
WRITE(IO)  DM,  GM 
ENTHY  EIELD 
REWIND  10 
READ(IO)  UM.  GM 
IE  (  ISYM.EU.O)  GO  TO  202 
DO  203  II  =  1.  LCMAXH 
203  VEL (LCMAX  ♦  1  -  II)  =  ISYM  *  VEL ( 1 1 ) 

202  CONTINUE 

IE  (NOPR.EQ. 1 )  GO  TO  1111 
88  IE  (WIN.NE.O. )  PRINT  744 
IE  (RIN.EQ.O. )  PRINT  746 
746  FORMAT ( 1 H 1  60X  1 9HC YL l NOER  TRANSDUCER//) 

744  FORMAT  (lHl  bOXlSHRING  TRANSDUCER//) 

PRINT  701  ,R  IN, ROUT, H.EK, NODI, I MAX, JMAX 
701  FORMAT ( *  RIN=«E7.3*  ROUT=*E7.3*  H=*F7.3*  FK=*E7.3»  N0ul=«I3 

l*  | mAX  =  * I  3*  JMAX  =  »13/) 

PRINT  36S 

36S  FORMAT  ( 9X JHNUM8X 7HSP  REAL11X12HSP  IMAGINARY9X8HVEL  REAL  10X1 3HVE 

SL  IMAGINARY/) 

1111  CONTINUE 

IF  (  ISYM.EO.-l )  GO  TO  99 
DO  7  J=l, LCMAXH 

VtLl(J)  =  .6* ( VtL ( J)  ♦  VEL (LCMAX  *  1  -  J)  ) 
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7  CONTINUt 

IF(MXU.NE.l)  GO  TO  806 
00  80S  J=l. LCMAXH 
IF (JMXT(J) .NE.O  )  GO  TO  80b 
00  807  I=1,LCMAXH 
DM ( I » J)  =  -  GM(I,JI 

807  GM ( i,J)  =  (0.0. 0.0) 

805  CONTINUt 

806  CONTINUE 

DO  246  1=1 .LCMAXH 
Vd)  =  (0.0. 0.0) 

DO  246  J=l. LCMAXH 

V(I)  =  V ( 1 >  ♦  GM ( I . J) *VEL 1 ( J) 

246  CONTINUt 

CALL  SIMXIDM. LCMAXH, V. ANSI  ) 

IF  (  ISYM.EO.  1  )  GO  TO  >78 
99  CONTINUt 

DO  IS  1=1. LCMAXH 
DO  IS  J=l. LCMAXH 

OM(I.J)  =  DM(I.J)  -  UMd.LCMAX  ♦  1  -  J)  - 
IS  GM(I.J)  =  GM(l.J)  -  GM ( 1 .LCMAX  ♦  I  -  J)  - 
DO  8  J=l, LCMAXH 

VELKJ)  =  .5*  (VEL  ( J)  -  VEL  (LCMAX  ♦  1  -  J) 

8  CONTINUE 

00  248  1=1, LCMAXH 
V(l>  =  (0.0, 0.0) 

00  248  J=l, LCMAXH 
V ( I >  =  V(I>  «  GM ( I , J ) * VEL 1 ( J) 

248  CONTINUE 

CALL  SIMX(OM, LCMAXH, V.ANS2) 

98  DO  10  J=l« LCMAXH 

ANS(J)  =  ANSl(J)  ♦  ANS2 ( J ) 

10  ANS (LCMAX  .  1  -  J)  s  ANSl(J)  -  ANS2(J> 

IF (MXD.NE. 1 )  GO  TO  809 
DO  808  J=l. LCMAXH 
IF (JMXT(J) .NE.O  I  GO  TO  808 
VEL ( J)  =  VEL (LCMAX* l-J)  =  ANS(J) 

ANS ( J)  =  ANS(LCMAX*1-J)  =  (0.0, 0.0) 

808  CONTINUE 

809  CONTINUE 

SUM J  =  (0.0, 0.0) 

DO  537  J= 1 , JM AXH 

IF (MIN. EG .0.)  GO  TO  600 

AREA ( J)  =  OELZ*R  IN 

AREA ( JMAXH, IMAX»J)  =  OELZ*ROUT 

GO  TO  S37 

600  AREA ( IMAX, J)  =  ROUT*OELZ 
#11  CONTINUE 

JRING  =  JM AXH 
IF (RIN.tO.O.)  JR ING=0 
DO  S38  J= 1 , IMAX 
AREA( JRING* J)=  HR ( J ) *UELH 
S3H  CUNTINUt 

00  539  J=l» LCMAXH 
S39  AREA (LCMAX  ♦  l  -  J)  =  AREA ( J) 

DO  7778  I  =  1,  LCMAXH 


UMd.LCMAX  ♦  1  -J) 
UMd.LCMAX  ♦  1  -J) 

) 
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7778  PURII)  =  ANS(I) 

00  401  J=1.LCMAX 

SUM  Ja  SUMJ  ♦  CON JG ( VEL  ( J)  )  *ANS I J) *ARE A  ( J) 

ANS  <  J  >  =  RHOC*ANS<J> 

IF(NOPR.EG.l)  GO  TO  1112 
PRINT  400  *  J* ANS  < J) «VEL  < J) 

400  FORMAT  C1X*I10«2C (E20.8.E20.8) 1 

1112  CONTINUE 

401  CONTINUE 

A1  =  RIN*ROUT 

AREAT  =A1*(  H  ♦  H  ♦  ROUT  -  RIN) 

SUMJ  =  -FKM0. 0*1.0)  *SUMJ/AREAT 
IFINOPR.EO.l)  GO  TO  1113 

PRINT  506.  SUMJ  _  „  „  .  . 

506  FORMAT!/*  THE  COMPLEX  RADIATION  IMPEDANCE  IN  UNITS  OF  RHU  C  A  IS  ( 
S*C(E15.8*E16.8)*)*) 

1113  CONTINUE 

402  CONTINUE 
END 
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SUROOUTTME  FR1KIG 

Type-  COMPLEX  GM.ISM.SIM.ITM.TIM.IIM.SSM.TSM.STM.SRM.TRM.DM.RHOC. 

lilisiiiliiilisiis- 

3  ,|Svn0.20».SIV(10.?0>,fTWI?0.1«».TIV(10.?»» 

.  iiuMft»20UPL7l?0)  »P7<?0) 

COMMON/RLKD/TSM/IIJ/T0M/P J/STM/BLKB/SSM/VEl  /  VS* V I/VELT/V ( 

COMMON/Pl T/JMAX. JMAXH. IMAX/RAn/RIN.R0UT/5/0M/6/GM. REST 
COMMON/TIDY/E ART/RCC/JTOP/RLKl/H, A, FK.P! 

FKA  =  FK* A  %  OFLZ  =  H/ JMAXH  *  OELR  =  (ROUT-R IN) / IMAX 

no  1  J=  1  «  JMA  XH 

KC  =  JMAX  ♦  1  -  J 

P2(j)  =  H  ♦  .S*OFLZ  *  J*0EL7 

PZ(KC)  =  “P7 ( J I 
PLZIM  =  P7  ( J)  -.S*OFL7 

1  PLZ  (KC 1  =PZ  (KC  1  -  .S*DEL7 
00  ?.  1  =  1  *IMAX 

2  pr ( I i  =  RIN  -  .S*OELR  ♦  I*OELR 
CALL  CALRES 

CALL  CALTRTG 
CALL  CTRM 
CALL  STTRIG 
CALL  CSTM 
CALL  SSTRIG 
CALL  CSSM 
CALL  TSTRIG 
CALL  CTRM 

00  *000  I=1*LCMAXH 
00  *000  J” 1 .LCMAXH 
GMU.J)  =  (0.0*0.01 
* 000  OM(I,J)  =  (0.0.0.01 

IF (RTN.FO.0. 1  GO  TO  RO 
00  *4  1  =  1  * IMAX 
00  *4  J=  l  .  IMAX 

GM (IMAXH* I  * JMAXH  *J)  =  TTV(I.J) 

*4  CONTINtJF 

00  31.  1=1. JMAXH 
KITE  =  JMAXH  .1-1 


00  3?  .1*1*  JMAXH 

OMU.Jt  =  IIM(KITF. JMAXH  .1  -  J) 
GM(I.J)  =  ItV(KlTE» JMAXH  .1  -  J) 

GM  LCMAX  .  1  -  J)  =  TIV(KITE, JMAXH 

32  OMniLCMAX  .  1  -  J.  =  I IM  (K I TE  .JMAXH 
no  33  1=1. IMAX 

GM(I.JMAXH.J)  =  ITV (K1TE.J) 

33  OM ( T .JMAXH. ))  =  ITM(KITE.J1 


no  34  1=1. JMAX 

GM(T. JMAXH. IMAX«J)=  ISV(KITF.J) 
34  OM ( I , JM A  XH* I M A  X ♦ J 1  =  ISM(KITF.J) 


♦ 

♦ 


J) 

J) 
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no  35  J=1  •  I max 

GM(I.JMAXH*LCMAXH*J>=  ITVMMAX  ♦  I  -K  I  TE  ,  I M  AX  ♦  1  -  J> 

35  DM  ( I  .  JMAXH*I_CMAXH*  J>  =  !TM<)MAX  ♦  1  -K  I  TE  .  I  MAX  ♦  1  -J) 
31  CONTTNUF 

DO  36  1=1 , IMAX 

KOOT  =  IMAXH* 1 

DM (KOOT.KOOT )  =  <0.5.0.01 

DO  37  .1=1  IMAXH 

K1TF  =  JMAXH  ♦  1  -J 

GM (KOOT «  J)  =  T I V  1 1 »K ITF ) 

DM (KOOT  »  J)  =  TIM(I.KITF) 

GM (KOOT.LCMAX  ♦l-J)=  T I V ( 1 . JM AXH ♦ J ) 

37  DM (KOOT.LCMAX  «1-J>=  T !M ( 1 . JMAXH* J) 

DO  30  J=1.JMAX 

GM (KOOT  , JMAXH*IMAX*J)=  TSV(l.J) 

30  OMIKOOT*  JMAXH*IMAX*J)=  TSM(l.J) 

DO  39  J= l » IMAX 

GM (KOOT.LCMAXH* JMAXH* J) =  TRV ( I » IMAX ♦ 1- J) 

39  DM (KOOT .LCMA XH* JMAXH* J >  =  TRM ( 1 ♦ ! MAX ♦ 1 - J* 

36  CONTTNUF 

00  40  1=1. JMAXH 

KOOT  =  JMAXH  ♦  IMAX  ♦  I 

00  41  J=l. JMAXH 

KIT  =  JMAXH* 1— J 

GM (KOOT. J)  =  STV(I.KITF) 

OM(KOOT.J)  =  SIMM, KITE) 

DM (KOOT.LCMAX* 1-J)=  SIM ( I « JMAXH* J) 

GM (KOOT.LCMAX* 1 - J »  =  SIV ( 1 .JMAXH* J) 

41  CONTTNUF 

DO  4?  J= 1 .IMAX 

DM(K00T.JMAXH*LCMAXH*J)=STM(JMAX*1-1.1MAX*1-J) 
GM (KOOT, JMAXH*LCMAXH*J) =STV<  IMAX* l -I . IMAX* 1“J> 
GM  (KOOT  ,JMAXH*J)=  STVU.J) 

42  DM (KOOT .JMAXH*J)=  STM(I.J) 

DO  43  J= 1  *  JM AX 

GM(KO0T.JMAXH*IMAX*J)  =  SSVU.J) 

43  DM (KOOT .JMAXH* IMAX*J)  =  SSM(T*J) 

40  CONTTNUF 
RFTURN 

00  DO  04  1=1 » IMAX 
00  04  J=  l.IMAX 
04  GM ( T . J I  =  TTV(I.J) 

DO  71  1=1, IMAX 
DM ( I  ,  I ) =  (0.5,0.01 
DO  7?  .1=  1 ,  JM  AX 
GM(T,IMAX*J)  =  TSV(T.J) 

72  DM ( T . IMAX* J)  =  TSM(T.J) 

DO  71  1=1. IMAX 

GM  ( I « JMAX  ♦  IMAX  ♦  J)  =  THVd.IMAX  ♦  1  -  J) 

71  DM ( I , JMAX  ♦  IMAX  *  J)  =  THM ( I  ,  IMAX  ♦  1  —  J > 

DO  73  1=1, JMAXH 
K  =  IMAX  ♦  1 
DO  74  J= 1 » IMAX 
GM  ( K  ,  J 1  =  STVU.J) 

74  OM(K.J)  =  STM(T.J) 

DO  75  J=1 .JMAX 
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GMIK.IMAX  ♦  J)  =  SSV(1.J> 

75  DM  (K «  IMAX  ♦  J)  =  SSMUtJ) 

00  71  J=1.IMAX 

GM ( K  »  JMAX ♦ !M AX ♦ J ) -  ♦  STViJMAX  ♦1-ItIMAX  ♦  1  -J> 
73  nM(K»JMAX*lMAX*J)=  ♦  5TM(JMAX  ♦l-JtIMAX  ♦  1  -J) 
ENO 
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SUBROUTINE  CTBM 
C0MM0N/FPS1 /EPS 

COMMON/S/BTO (10, 10,32)  ,BTU  1 1  , 10,32) »QZ  (32, 10) .SIGMA (32) .REST 1128) 
COMMON/G/GTO (10.10,3?)  ,GT1  (11 .10,32) .RUST1  (2880) 
common/rlki/h.a.fk.pi/pit/jmax.jmaxh.imax/ttv/ttv/tidy/east 
COMMON/RLK2/X I » WF/TOY/NQD1  .N2.N3/BLK3B/BI O.RI 1 /BLK3/BS0  »BS1 

COMMON/RLK A/PR.PLR.PZ.PLZ/OEL/DFLR.OELZ/BOLD/GSO  »GS1 /RCC/ JTOP 

COMMON/FS/FKFK ,FKF  AST  .  TEK/TBV/TBV/UJ/TRM/GGMS/GGMS 

DIMENSION  P7 (20 ) ,PL7 (20 ) ,PR ( 20 ) ,PLR ( 20 ) , TTV ( 10 , 1 0 ) »TBM ( 1 0, 1 0 ) , 

1  TBV(IO.IO) .GGMS (10.32) 

TYPE  COMPLEX  TTV.TRM.TRV 

DIMENSION  XI (32) • WE (32) ,SUMT (10) ,GS0 ( 10.32)  ,GS1 ( 1 0 .32) ,UMT ( 1 0 ) » 
IRUMT(IO),  RSOOO. 32) .RSI (10,32) .BIO (10,32) .RI1 (10,32) 
common/storsctb/tbsin.trcos.tbexp 
OIMFNSION  TBSIN(32) ,TBCOS(32) .TBEXP (32, 10) 

EQUIVALENCE  (TBCOS ( 1 ) .TBEXP ( 1  )  ) 

TPK  a  0.2S  «  FKF AST 

PIOTWO  =  O.S  •  PI  S  TWOPI  =  PI  ♦  PI 
N=N001 

HO  a  ,5*0ELR 
ARGP  a  H  ♦  H 
00  7  JRal.lMAX 
R2  =  PR (JR) 

RL  =  R2  -  HO 
RU  a  R2  ♦  HO 
DO  7  IRal.lMAX 
R 1  a  PR(IR) 

UMIaSUMR=SUMI=SUMRB=SUMIB=0. 

00  1  L  =1«N 

YAK  a  RTO(IR.l.L)  •  BTl(JR.l.L) 

VI  a  YAK  *  TBCOS (L  > 

VR  =  YAK  *  TBSIN(L) 

SUMTB  a  SUMIB  -  VI 
SUMRR  a  SUMRB  ♦  VR 
SIIMR  a  SUMR  -  VI  •  SIGMA(L) 

SUM l  a  SUM)  ♦  VR  *  SIGMA (L ) 

1  IJMI  a  UMI  ♦  YAK 
UMI  a  UMI»TFK 

R  =  PI 

IF (R2.LE .R1 )  R  =  R2 
IF(R.EO.O.)  R  =  OELR 
FACT  =  2./(R*PI*ARGP> 

RM  a  PI 

IF(P.EO.Rl)  RM  a  R2 

FACT  1  a  PIOTWO  /  ABS(ABS(R?-R1)  -  HO) 

EFACT1  a  1.0  -  FXPF(-FACT1  *  ARGP) 

UMT(I)  a  0.  %  SUMT(l)  a  -SUMP  $  RUMT ( 1 )  =  -  SUMRR 
FL  =  0.0 
ICE  =  1 

2  ICE  a  ICE  ♦  1 

UMMTT  a  SUMMIT  =  SUMMITS  =  0. 

DO  37S  Ma 1 ,N 

YAKR  =  BTO(IR.ICE.M)  •  BT1 ( JR, ICF.M) 

VAKR  a  YAKR  -  GTO(IR.ICE.M)  *  GT 1 ( JR » ICE ,M) 

UMMTT  =  UMMIT  ♦  VAKR 


00009400 


00011300 
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ART T STB  =  YAKR  *  TBEXP(M»ICE) 

ARTIST  =  ARTISTB  *  QZ(M«ICE) 

SUMMIT  =  SUMMIT  ♦  ARTIST 
175  SUMMITB  =  SUMMITB  ♦  ART ISTB 

SUMT(ICE)  =  SUMT(ICE-l)  ♦  SUMMIT 
BUMT(ICF)  =  BUMT(ICF-l)  ♦  SIJMMITB 
UMT (ICE)  =  UMT(ICE-l)  ♦  TPK  *  UMMIT 
FL  =  FL  ♦  FKFAST 

FRFC  =  ABSIFACT  *  RM  *  EXPFI-ARGP  *  FL )  *  EFACTI  /(FL*SUMT (ICE) )  ) 
TESTFR  =  ABS (1.0  -  SUMT (ICE) /SUMT ( ICE-1  )> 
if (tester.lt .eps.ano.erfc.lt .EPS)  GO  TO  3 

IF(TCE.GE.JTOP)  GO  TO  3 
GO  TO  2 
3  CONTINUF 

ANSR  =  SUMT (ICE)  %  ANSI  =  -SUMI 

VNSI  =  -SUMIB  «  VNSR  =  -  BUMT(ICE) 

TBM ( IR» JR)  =  CMPLX (ANSR* ANSI )  S  THV(IP.JR)  s  CMPLX (VNSR t VNSI ) 

OK  =  4 . *R 1 *RU/ ( ( R 1  ♦  RU)  *  (R1  ♦  RU) ) 

CU  =  ( R 1 *RU) *ELL IPE (DK )  ♦  (PII-RI ) *ELL IPK (DK ) 

DK  =  4.*R1*RL/((P1  ♦  RL)  *  (R1  ♦  RL)) 

CL  =  (R) *RL) *ELL IPE (DK)  ♦  (RL-R1 ) *ELL IPK (OK) 

CF  =  (CU  -  CL)  /  TWOPI 
ANSR  =  -  UMT (ICE)  -  CF 
7  TTV ( TR» JR)  =  CMPLX (ANSRfUMI ) 

ENO 
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SURPOUTTNE  CSTM 

COMMON/ICE/ICE.NOFXJ 

TYPE  COMPLEX  STV.ITV.STM.ITM 

COMMON/EPS/EPS/ JOE/ JTOP 

COMMON/STV/STV/1TV/ITV/PJ/STM/ITM/ITM 

COMMON/RAD/RIN.ROUT 

DIMENSION  ITM ( 20 . 1 0 ) .  I  TV (20 . 1 0 ) .STM (20 . 1 0 > ,STV (20, 1 0 > 

COMMON/PLK A/PR .PER, PZ.PLZ/PIT/JMAX.JMAXH.IMAX 
DIMENSION  PZ(20)  ,PL7(20> ,PR<20> ,PLR(20> ,GGMS(10,32) 

DIMENSION  RSOUO.3?)  ,RS1  ( 1 0 , 32  >  »RI0  ( 1 0 . 32)  ,RI1  (10,32) 

DIMENSION  XI (32), WF (32), SUMV (10) 

DIMENSION  GSO (1 0,3?)  ,GSI (10,32) .RUMT ( 1 0 ) ,SUMT < 1 0 ) ,BUMV ( 1 0 > 

TYPE  COMPLEX  ANS«ANSB»VNS»VNSB 
COMMON/FS/EKFK ,EKFAST, TFK/GGMS/GGMS 

COMMON /DEL /DELR, DEL 7/ROL D/GS 0 ,GS1 /RLK 1 /H, A,EK,PI/Bt.K?/X I ,WE 
COMMON/S/BTO ( 10,10,3?)  ,BT1  (1) ,10,32) ,07(32,10) .SIGMA (32)  .REST (128) 
COMMON/TOY/N»NN»NNN/BLK3B/BIO,B!1/BLK3/BSO,BS1/T!DY/FAST 
COMMON/B/GTO (10,10.32) .STEXP ( 32,20, 1 0 ) 

COMMON/STORSCST/STSIN  ( 32,  ?1 ) .STCOS ( 32,21 ) 

PIOTWO  =  O.S  •  PI 
PI  =  A 

HD  =  „5*0ELR 

DO  1000  JR= 1,1 MAX 
R2  =  PR ( JR) 

RL  =  R?  -  HO 

RU  =  R?  ♦  HD 

DO  1000  !R=1,JMAX 
ARGP  =  H  -  P7(IR> 

SUMR=Sl)MI=SUMRR=SUMIB  =  VUMR=VUM I=VUMRB=VUM IR=  0. 

DO  1  L  = 1 «N  00009400 

BU  =  BT 1 ( JR ,  1  ,L  ) 

VAK  a  RSOO.L)  *  BIJ 
VAKR  =  RIO(l.L)  *  BU 
YAK  =  VAK  *  SIGMA (L ) 

YAKR  =  VAKR  *  SIGMA  (L ) 

SUMP  =  SUMR  -  YAK  *  STCOS (L,IR) 

SUM!  =  SUMI  ♦  YAK  «  STSIN(L.IR) 

VUMR  =  VUMR  -  VAK  »  STSIN(L.IR) 

VUM I  =  VUMT  ♦  VAK  *  STCOS(L.IR) 

SUMRR  =  SUMPR  -  YAKR  *  STCOS (L, IP) 

SUM T R  =  SUMIR  ♦  YAKR  *  STSIN(L.IR) 

VUMPR  =  VUMRR  -  VAKB  *  STSIN(L.IR) 

VUM I R  =  VUM I R  ♦  VAKR  *  STCOS(L.IR) 

1  CONTINIJF 
R  =  O) 

IF(P?.LE.R1)  R  =  R2 
IE(P.EO.O.)  R  =  DELP 
FACT  =  ?./(R*PI*APGP) 

RM  =  PI 

IF (P.EO.R1 )  RM  =  R? 

FACT  1  =  PIOTWO  /  AMTN1 (POUT  -  R?,  R2  -  RIN) 

EE  ACT  1  =  1.0  -  FXPF (-APGP  *  F  ACT  1 ) 

SUMT(l)  =  -  SUMR 
RUMT  ( 1 )  =  -  SUMRR 

SUMV(l)  =  VUMR  S  RUMV(t)  =  VUMRB 
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Ft  =  0.0 

ICF  =  1 

2  ICF.  =  ICF  ♦  1 
VUMMIT  =  VUMMITR  =  0. 

SUMMIT  =  SIJMMITR  a  0. 

DO  175  M=1,N  00011300 

RU  =  STFXP(M.IR.ICF)  o  RTl  ( JR.  ICF. Ml 
V AKR  =  RSOUCE.M)  •  BU 
VAKRR  =  RIO(ICF.M)  *  Rl) 

SUMMIT  =  SUMMIT  ♦  VAKR  *  OZ(M.ICF) 

SUMMTTR  =  SUMMITS  ♦  VAKRR  *  OZ(M.ICF) 

VUMMIT  a  VUMMIT  ♦  VAKR 
VUMMTTB  =  VUMMITR  ♦  VAKRB 
375  CONTTNUF 

RUMV  ( ICF )  =  RUMV ( ICF  -  11  ♦  VUMMITR 

SUMV(ICF)  =  SUMV ( ICF  -  1)  ♦  VUMMIT 

SUMTUCF)  =  SUMTUCF  -  1)  ♦  SUMMIT 

RUMT(ICF)  =  BUMT ( ICF  -  1)  ♦  SUMMITS 

FL  =  FL  ♦  FKFAST 
FLT  =  FI  ♦  FACT  1 

ERFC  =  ARSIFACT  *  RM  *  FF ACT  1  *  FXPFl-ARGP  *  FL)  /  <FL*SUMT ( ICF) )) 

IF (FBFC.LT .FPS)  GO  TO  3 
IF(ICE.GE.JTOP)  GO  TO  3 
GO  TO  ? 

3  CONITNUF 
ANSI  =— SUM  I 
ANSR  =  SUMT(ICF) 

ANSIR  a  -  SUMIR  %  ANSRB  =  BUMT (ICF) 

NOFX I  a  0 

CALL  CF INT ( ARGP.  RU.  POUT,  FL.  CFU) 

CFST  a  RU  *  CFU 
NOFX I  a  1 

CALL  CFONEUBGP.  RU.  POUT.  FL .  CFU) 

CFSTV  a  RU  *  CFIJ 
NOFX I  a  0 

CALL  CFINTIARGP.  PL.  POUT.  FL.  CFL) 

CFST  a  CFST  -  (PL  *  CFL) 

NOFX I  a  I 

CALL  CFONE ( ARGP.  RL .  ROUT.  FL  .  CFL) 

CFSTV  a  CFSTV  -  (RL  *  CFL) 

NOFX I  a  0 

CALL  CFTNT ( ARGP.  PU.  PIN.  FL •  CFU) 

CF I T  a  RU  *  CFU 
NGFX 1  a  1 

CALI  CFONE ( ARGP .  RU.  PIN.  FL »  CFU) 

C F I T v  =  RU  *  CFU 
NOFX I  a  0 

CALt  CFTNT ( ARGP.  RL.  RIN.  FL ♦  CFL) 

CF I T  a  CFIT  -  (PI  *  CFL) 

NOFX I  a  1 

CALL  CFONF(ARGP.  »L.  RIN.  FL.  CFL) 

CF  I  TV  a  CFITV  -  (Rl.  *  CFL) 

ANSR  a  ANSR.CFST  $  ANSRB  a  AN5RR  ♦  CFIT 
STM(IR.JR)  a  CMP|.  X  (ANSR.  ANSI ) 

ITM(IR.JR)  =  CMPI X (ANSRR. ANSIR) 

SUMV(ICF)  a  SUMV ( ICF )  ♦  CFSTV 
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RUMV(ICF)  =  RIIMV(ICF)  ♦  CFTTV 

VMS  =  CMPL X (-SUMV ( ICF ) •  VUMI )  S  VNSB  =  CMPLX (-BIJMV ( ICF ) « VUM IB) 
STV ( TR* JR)  =  VNS  t  ITV(IR.JR)  =  VNSR 

inoo  conttnuf 
ENO 
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SURPOUTINF  CSSM 
COMMON/ T  CF/ICF.NOEX i 

0  IMF  MS  TON  STV(  10,20)  ,119(10,20)  .SSVH0.20)  ,  ISV(  10,20) 

OIMFNSION  S IM (1 0 ,20 ) .  I IM (1 0 ,20 >  . SSM (1 0 ,20 ) , I SM (1 0,20 > 

TYPF  COMPLFX  SSV.SIV.ISV.tlV.SSM.SIM.ISM.HM 
OIMFNSION  P7(20)  ,PL7(?0) ,PR(20) .PLM20) 

TYPF  REAL  ITIV.TIRV.ISIV.ISRV.ISR.IIR.IIT.IST.TII.IIR.ISI.ISR 
OIMFNSION  SSR(IO) ,  IIBdO)  .SIROO)  ,  TSRdO) 

OIMFNSION  GSO(10,3?) »GS1 ( 1 0 , 3? ) » X  I  (32f ,WF(32) .SUMT(IO) 

OIMFNSION  RS0(10,3?> .RSI (10,32) 

OIMFNSION  RIO (1 0,32) ,RI1 (10,32),GIO(10,32>,GI1 (10.32) 

TYPF  COMPLFX  SS.II.YAKC.UMII.UMSS.UMIS 

OIMFNSION  SST(10),I1T(10).SIT(10)  ,IST(10)  .GGMSU0.3?) 

OIMFNSION  RIIMT  ( 10)  .OIIMT(IO)  ,EXPF1  (32) 

COMMON/FPS/FPS 

COMMON/PAD/RIN,ROUT/TOY/NOOl.NQO?.NQD3/J0F/JTOP/R0LD/GS0,GSl 
COMMON/P  I T/JMAX, JMAXH,IMAX/RLKA/PR,PLR,PZ. PI  Z 
COMMON/SS V/SSV  / 1  IV/IIV/SIV/SIV/ISV/I SV/GGMS/GGMS 
COMMON/RLKR/SSM/SIM/SIM/I IM/I IM/ISM/ISM/FS/FKFK,FKFAST,TFK 
COMMON/RLKl/H,A,FK.PI /RLK2/XI .WF/BLK3/RS0 . RS 1 /T IOY/F AST 
COMMON /RROLO/G I0,GI1/BLK3R/B10,RI1  /OFL/OELR  »DELZ 
COMMON/5/RTO(10.10,32),RT1  (11,10.32) »QZ(32»10> .SIGMA (32) .REST ( 128) 
COMMON/(S/GTO(  10. 10.32)  .SSEXP ( 32,20 . 1 0 ) 

COMMON/STORSCST/SSS  IN ( 32.21 ) ,SSC0S(32,21 ) 

N  =  NOOl 

APGP  =  ,S*OFLZ 

7ERO  =  0.0 

TPK  =  . ?5*FKF  AST 

70  =  H-  ,5*OEL7 

SSI=SSR=I 11  =  1 IR=SI  I  =  SIR=ISI  =  ISR=0. 

IJMTI  =  I  IMIS  =  IIMSS  =  (0,0. 0.0) 

00  1  L  =1 »N  00009400 

SIGMAP  =  ARGP*S IGMA (L ) 

EXPFl(L)  =  FXPF(-07(L,1)*ARGP) 

TRIGFR  =  2. *COS (SIGMAP)  -  2. 

TRIGFI  =  2 . *S IN (S I GMAP) 

YAK  =  TFK»WF (L)/SIGMA(L) 

YAKP  =-YAK*TRIGFR 
YAK  I  a  YAK*TRIGFI 
YAKC  =  CMPLX (YAKP.-YAKI) 

SSI  =  SSI  ♦  R  ,0<1,L)*RS1 (1«L)*YAKI 
SSR  =  SSR  ♦  RS0(1,L)*RS1(1.L)*YAKR 
HR  =  MR  -  RI0(1,L)*BI1  (l.L)*YAKR 
III  a  HI  -  RI0(1,L)*RI1(1.L)*YAKI 
SII  a  SII  -  RS0(1.L)*RI1(1.L)*YAKI 
SIR  a  SIR  -  RS0(l,L)*R!l (1,L)*YAKR 
ISI  a  ISI  ♦  BI0(1,L>*BS1(1,L>*YAKI 
ISR  a  ISR  ♦  B 1 0 ( 1 »L ) *RS1 (1,L)*YAKW 
UMTI  a  UM 1 1  ♦  YAKC*BI0(1.L)*BI0(1.L) 

UMTS  =  UM I  YAKC*B10(1,L)*RS0(1 .L) 

IJMSS  =  UMSS  ♦  YAKC*RS0(1 .L)«RS0 ( 1 ,L> 

1  CONTINUE  00010900 

SUMTd)  a  REAL  (UMSS)  %  RUMTO)  a  REAL  (UMI I )  %  OUMT  ( 1 )  a  REAL  ( UM  IS) 

SST  ( )  )  a  SSP  1i  I  I T  (  1  )  a  HR  %  SIT(l)  a  SIP  %  1ST  (  1  )  a  ISR 
FXF ACT  a  EXPF (-FKF AST*ARGP ) 
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no  sm  ICE  =  2,  1  i 
SUMSS=SUMII=SUMSI  =  SUMIS  =  n. 

SUMMIT  S  SUMMIT^  =  SUMMITD  =  0. 
no  175  M=1,N 

FXPFl(M)  s  FXPF1  (M)*FXFACT 
EXPF?  =  2.* (EXPF ] (M)  -  1.) 

YAK  =  -TPK*WF(M> 

YAKP  =  -EXPF2*YAK 
YAKRG  =  YAKP  /  OZIM.ICF) 

YAK  =  YAK  ♦  YAK 

SUMSS  =  SUMSS  ♦  RSO  I  ICF,M)  *  RS1IICF*M)  *  YAKRG 
SUMSS  =  SUMSS  -  GS0(ICF,M)*GS1  (ICE,M)*YAKR 
SIJM 1 1  =  SUM  I  I  -  RI0IICF,M)  »  HIl(ICFfM)  *  Y  AKRG 
SUMIT  =  SUM  1 1  ♦  GI0(ICE»M)*GI1  (ICE»M)*YAKR 
SUMS!  =  SUMSI  -  RSO(ICF.M)  »  Bll(ICE*M>  *  YAKRG 
SUMS!  =  SUMSI  ♦  GS0(ICF,M)*GI1  (TCE,M)*YAK 
SUMIS  =  SUM  VS  ♦  R 1 0  I  ICE »M )  *  RSI ( ICF  *M)  *  YAKRG 
SUMIS  =  SUMIS  -  GI0(ICF,M)*GS1 (ICE,M)*YAK 
RSRS  =  RSO(ICF.M)*RSO(ICF,M) 

RIRS  =  HIOIICF  iM)*RSO(ICF»M) 

RIRT  =  RIO(ICE*M)*BIO(ICE»M) 

SUMMIT  =  SUMMIT  ♦  RSRS  *  YAKRG 
SUMMITS  =  SUMMITR  ♦  HIRI  *  YAKRG 
SUMMTTO  =  SUMMITR  ♦  SIRS  •  YAKRG 
775  CONTINUE 

SSTITCE)  =  SST< ICE-1)  ♦  SUMSS 
IIT(ICF)  =  IlTIICE-l)  ♦  SUM  1 1 
SIT ( ICF)  =  SIT(ICF-l)  ♦  SUMSI 
1ST ( ICF)  =  1ST  <  ICF  -1>  ♦  SUMIS 
SUMT(ICF)  =  SIJMTUCF  -  1)  ♦  SUMMIT 

RUMT I ICF)  =  RUMT  ( ICF  -  1)  ♦  SUMMITR 
SOI  OUMT(ICF)  =  OUMTITCF-I)  ♦  SUMMITR 
AA  =  4.*ROUT*ROUT 
P  =  ARGP 
PP  =  P*P 

OK  =  AA/IPP  ♦  A A ) 

SDK  =  SORT (OK) 

CF2  =  -P*SOK*ELLIPK(DK)/(PI*AA> 

SSR=  ROUT* (SST  I  JTOPt  ♦  CF2  ♦  CF?) 

AA  =  RIN*RIN*4. 

OK  =  AA/IPP  ♦  A A ) 

SDK  =  SORT (OK) 

CF?  =  -P*SDK*FLLIPK(DK)/(PI*AA) 

IIP=  RIN* 1 1  IT  I JTOP)  -  CF?  -  CF?) 

SSI  =  ROUT*SSI  %  1 1 1=  R  IN*  1 1 1  S  SII=  RIN*SI I  t  ISI*ROUT*ISI 
ICE  =  10 

FL  =  (ICE-1 >*FKFAST 
NOEXT  =  0 

CALL  CF INTI ARGP. ROUT, RIN, FL.CF IS) 

CALL  CFINTIARGP. RIN, ROUT, FL,CFP)  $  CFS1=  -CFP 

SIT  I ICF)  =  SIT  I  ICE)  ♦  CFSI  ♦  CFSI 

1ST  I ICF)  =  1ST (ICF)  ♦  CFIS  ♦  CFIS 

SIR  =  RIN*SIT  I  ICF)  %  ISR  =  ROUTMST  I  ICE) 

1  SR  =  ISR  -  1.0 
JM  =  <3 

FLAM  *  FKFAST*JM 
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00013600 

00013500 
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CALI  CF7FR0I ARGP,ROUT»R0UT.FI  AM,CP) 

CALL  CFZFROIJ FRO, ROUT .ROUT. FLAM. CM) 

CF1  =  2MCP-CM) 

CALL  CF7FRO( ARGP.RIN.BIN, FLAM, CP) 

CALL  CF7FRO(7FRO,RIN,RIN,FLAM,CM) 

CF2  =  ?*(CP-CM) 

CALL  CF7FRO(ARGP.ROIIT.RIN.FLAM,CP) 

CALL  CF7FRO (ZFRO. ROUT.  R1N, FLAM, CM) 

CF3  =  -2MCP-CM) 

ANSR  =  SUMT(ICF)  ♦  CF1 
ANSRP  =  BUMT(ICF)  ♦  CF? 

ANSRO  =  IKJMT ( TCF)  -  CF3 
01  =-A IMAG (  UMII) 

IIV(l.l)  =  R1N*CMPLX  (ANSRB.OH 
01  =-AIMAG(  IJMIS) 

ISVH.l)  =  ROUT*CMPLX (ANSRO, 0! ) 

SIV(1,1)  =  RIN*CMPLX (ANSRO, Q!  ) 

01  =-A IMAG (  UMSS) 

ssvn.n  =  rout*cmplx(ansr,ot> 

IIM(1,1)  =-CMPLX(II^,I!T)  ♦  (0.5,0.01 
SSM  (1,1)  =  -CMPI.X  (SSR,SS!)  ♦  (0.5, 0.0) 

SIM ( 1 , 1 )  =  -  CMPLX  (SIR, SID 
I SM (1.1)  =  -CMPL  X(I5R,ISI) 

00  65  J=2, JMAX 
7  =  PL7(J) 

ARGP  =  70  -  Z  -  0EL7 
ARGM  =  70  -  Z 

SSI=SSR=III=IIR=SII=SIP=ISI=TSR=0. 
SSTV=SSRV=ITIV=IIRV=SIIV=SIRV=ISIV=ISRV=  0. 

00  61  L=1«N 

YAKR  =  -SSCOS  <L , J-l ) 

YAKI  =  SSSIN(L.J-l) 

I  SR  =  TSR  ♦  RIO(l,L)*RSm,L)*YAKR 

ISI  =  ISI  ♦  BI0(1,L)*BS1 ( 1 »L) *YAKI 

SIR  =  SIR  -  RS0(1,L)*RI1(1,L)*YAKR 

SII  =  SII  -  RS0(1,L)*R11 (1 ,L)*YAKI 

III  =  III  -  BI0(1,L)*RI1(1«L)*YAKI 

HR  =  IIR  -  RI0(1»L)*BI1U,L>*YAKR 

SSR  =  SSR  ♦  RS0( 1 ,L)*RS1 (1,L)*YAKR 

SSI  =  SSI  ♦  BS0(1,L)*RS1(1.L)*YAKI 

1 1 IV  =  IIIV  ♦  RTO(t'L)  *  BTO(l.L)  *  YAKI 

I IRV  =  IIRV  ♦  BIO(l.L)  *  RIO(l.L)  •  YAKR 

I SRV  =  ISRV  ♦  B 1 0 ( 1  ,L )  *  RSO(l.L)  *  YAKR 

ISIV  =  ISIV  ♦  BIO ( 1 ,L)  •  BSO(l.L)  *  YAKI 

SSRV  =  SSRV  ♦  RSO(l.L)  «  BS0(1,L)  *  YAKR 

61  SSIV  =  SSIV  ♦  PSO(l.L)  »  BSO(l.L)  •  YAKI 

SI  IV  =  ISIV  «  SIRV  =  »SRV 
F ACT  1  =  PI  *  0.5  /  (ROUT  -  RIN) 

FACT  =  ? . / ( A*P I  * ARGM ) 

SST(l)  =  SSR  %  I T T ( 1 ) =  IIR  5  SIT(l)  =  SIR  S  1ST ( | ) =  ISR 
SSR ( ) ) =  SSRV  %  I  IB ( 1 ) =  IIRV  %  SIB(l)  =  SIRV  %  TSB(l)  =  ISRV 
ICF  =  1 

2  ICF  =  TCF  ,1 

SUMSS  =  SUM  I  T  =  SIIMS  I  =  SUM  IS  =  0. 

VUMSS  =  VUM 1 1 =VUMS I  =  VUM I S  =  0. 
no  376  M=1,N 
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VAKR  =  SSFXP <M,J-I . TCF) 

VUMSS  =  VUMSS  ♦  RSO < 1CF.M)«RS0 ( !CE,M)*VAKP 
VUMIS  =  VI IM  I S  ♦  BIO(TCF,M)«RS0(!CE.M)*VAKR 
VIJMII  =  VIJMT!  ♦  BJO(  tCF,M)«B!0(1CF.,M)*VAKR 
SUMS  I  =  SUMS!  -  RS0(1CF,M)  *  R!1(ICF,M)  *  VAKR 

SUMfS  =  SUM  IS  ♦  R  I  0  <  I CF  *  M )  #  BSHICF.M)  *  VAKR 

SUMSS  =  SUMSS  ♦  RSO(ICF.M)  «  BSHICF.M)  *  VAKR 

376  SUM  1 1  =  SUM  I  I  -  RIOUCF.M)  «  R  M  (  ICF  «M)  *  VAKR 


VUMST  =  VUMTS 

SST(ICF)  =  SST(ICF-l)  ♦  SUMSS 
1 1 T  ( ICF  )  =  TITHCF-l  >  ♦  SUM!! 

1ST ( ICF )  =  1ST (ICE  -1)  ♦  SUMIS 
SIT ( ICF)  =  SIT(ICE-l)  ♦  SUMS! 

SSR(ICE)  =  SSRUCF-1)  ♦  VUMSS 
1 1 R ( ICF )  =  I !B ( ICF-1 )  ♦  VUM 1 1 
ISR(ICE)  =  ISR ( ICF  -1)  ♦  VUM  I S 
SIR(ICE)  =  SIR! ICE- 1 )  ♦  VUMSI 
FI.  =  (ICE  -  1  )  *FKF  AST 
FIT  =  FI  ♦  F  ACT  1 

FRF  C  =  (FACT/FL > « (FXPF <-ARGP*FL >  -  FXPFI-  ARGP*FLT ) ) 

FRFC  =  FRFC/S I T ( ICF ) 

FRF  C  =  ARS(FRFC) 

IF(FRFC.LE.EPS)  GO  TO  3 
IFUCF.GE.10>  00  TO  3 
00  TO  ? 

3  SSI  =  ROUT*SS  I  f.  I  1 1  =  R IN* ! I !  *  Sl!=  RIM*SII  5  ISI=ROUT*ISI 
NOEXT  =  0 

CAt.l.  CFINT (ARGM.  PIN.  ROUT.  FL .  CFM> 

CALL  CFTNHaRGP.  RIN.  POUT.  FL.  CFP> 

SITUCF)  =  SITUCE)  -  (CFP  -  CFM> 

CALL  CF!NT(ARGP.  RIN,  RIN.  FL,  CFP> 


NOFXT  =  1 

CALL  CF7ER0 I ARGP ,  RIN,  PIN,  FL,  CP) 
NOFXT  =  0 

CALL  CFINT (APGM,  RIN,  PIN.  FL,  CFM) 
NOFXI  =  1 

CALL  CF7FR0 ( ARGM ,  PIN,  PIN.  FL,  CM) 

1 1 T ( ICF )  =  ITT(ICE)  -  (CFP  -  CFM) 
IIBUCF)  =  IIR(ICE)  ♦  (CP  -  CM) 

NOFXI  =  0 

CALL  CFTNT (ARGP,  ROUT.  ROUT,  FL,  CFP) 
NOFXT  =  1 

CALL  CF7EP0(APC,P,  POUT,  POUT,  FL ,  CP) 
NOEXT  =  0 

CALL  CFINT ( ARGM,  ROUT,  POUT,  FL »  CFM) 
NOFXT  =  I 

CALL  CF7FP0(ARGM,  ROUT,  ROUT,  FL,  CM) 
SST(TCF)  =  SST(ICF)  ♦  (CFP  -  CFM) 
SSB(TCF)  =  SSfl(ICF)  »  (CP  -  CM) 

NOFXI  =  0 

CALL  CFINT (ARGP,  POUT,  RIN,  FL,  CFP) 
NOFXT  =  1 

CALL  CF7ER0 ( ARGP ,  ROUT,  RIN,  FL,  CP) 

NOFXI  s  0 

CALL  CF I  NT ( ARGM ,  POUT  ,  PIN,  FL,  CFM) 
NOFXI  =  1 
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CALL  CF7fC0( ARGM,  ROUT.  RIN,  FL  .  CM) 

IST(ICE)  =  IST(ICF)  ♦  (CFP  -  CFM) 
lSH(TCF)  =  ISB(ICE)  ♦  (CP  -  CM) 

SIR(ICF)  =  SIBOCE)  ♦  (CP  -  CM) 

SSR  =  ROUT*SST(ICE>  5  UR  =  RIN*IIT(ICE> 

SIR  =  RIN*S I T ( ICE )  %  ISR  =  ROUT*IST < ICE) 

SSIV  =  ROUT*SS  I V  S  I  I  IV  =  RINMIIV  *  SIIV  =  R1N*SIIV 
ISIV  =  ROUT  *  IS  IV 

SIRV  =  RIN*S  IB ( ICE )  $  ISRV  =  ROUT* ISB ( ICE ) 

SSRV  =  ROUT*SSB (ICE)  t  IIRV=  RIN*IIB(ICE) 

SS  =  CMPLX (SSR.SSI )  S  II  =  CMPLX ( I IR, 1 1  J  ) 

SSMO  «.I)=SS 

ISM ( I  * J) =  CMPLX ( I  SR* IS  I ) 

SIM(I»J)=  CMPLX  (SIR.ST I ) 

IIMO  .  J)=II 

SSV ( 1  * J) =  -CMPLX  (SSRV. SSIV) 

IIV(1.J)=  -CMPLX  (IIKV.ITIV) 

ISV(1.J)=  -CMPLX  ( ISRV*  ISIV) 

SIV(1.J)=  -CMPLX (SIRV. SIIV) 

IF ( J.GT . JMAXH)  GO  TO  66 

5SM ( J» I ) =SSM ( 1  * J)  5  SSV(J.1)=SSV(1.J>  S  ISM(J*1)=ISM(1»J) 

simu,i)=sim(i,j)  $  isv(j.n=isvu.j)  s  iiM(j.n  =  iiM(i ,j» 
iiv(j.i)  =  nvn»j)  %  siv(j,n=sivd,j) 

66  CONTINUF 

00  I ?  J= 1 . JMAX 

JILT  =  JMAX-J.l 

00  U  M=2» JILT 

IF (M.GT. JMAXH)  GO  TO  6G1 

SSM(M.M.J-l)  =  SSM(I.J) 

SSV (M. M ♦ J— 1 )  =  SSV ( 1 » J) 

I SM ( M  *  M*  J- 1 )  =  ISM(l.J) 

ISV(M,M»J-1)  =  ISV(l.J) 

IIM(M,M*J-1)  =  IIM(l.J) 

I IV (M.M* J-1 )  =  IIV(l.J) 

SIM(M,M*J-1)  =  SIM(l.J) 

SI  V (M.M* J-1 )  =  SIV(l.J) 

601  CONTINUF 

IF(M*J-1.GT. JMAXH)  GO  TO  11 
SSM(M»J-1,M)  =  SSM(l.J) 

SSV (M*  J— 1 .M )  =  SSV(l.J) 

ISM(M.J-l.M)  =  ISM(l.J) 

I S V ( M ♦ J— l.M)  =  ISV(l.J) 

IIM(M*J— l.M)  =  IIM(l.J) 

IIV(M.J-l.M)  =  IIV(l.J) 

SIM ( M*  J- l.M)  =  SIM(1 ,J) 

SIV(M.J-l.M)  =  SlV(l.J) 

11  CONTINUF 

12  CONTINUF 
FNO 
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subroutine  ctsm 

C0M''oN/l  CE/  ICE  i  NOE  X I 

TYPE  COMPLEX  T I  V . TSV . T TM. TSM  .  ANS. ANSR . VMS . VNSR 
DIMENSION  TSM(10.20)*TIM(lO»?0).TSV(10»?0)»TIV(10»20) 

DIMENSION  PZ(21t  *PL7<20) .PR < 20 ) .PLR < 20 ) 

COMMON/EPS/EPS 'jOE/JTOP 

COMMON /BLK ] /H.rt.FK.PI /RAD/RIN.ROUT /PIT/JMAX  « JMAXH* IMAX/IJK/I.J 

COMMON/TIM/TIM/RLKO/TSM/TIV/TIV/TSV/TSV/RLKA/PR.PLR.P7.PLZ 

DIMENSION  BSO Cl  0.32) .RSI (10, 321. EX (321 

DIMENSION  BIO (10.321.BI1  1 1 0*32). G 10 <10*321  ,GI1 (10.32) 

DIMENSION  SUMTV(IO)  .RUMTVC  1) .XI (32)  »WE  (321.SUMT  CIO)  *BUMT  (101 
DIMENSION  GSO (10*32) .GS1 ( 1 0 , 32) ,GGMS ( 1 0. 32) 

COMMON /S/BT  0(10*10*32) .BT  1  ( 1 1 «  10.32) ,QZ  <32*  1  0) .SIGMA (32) .REST ( 1 
COMMON/6/Gf 0(10.10. 32) .TSE XP ( 32. 20 . 1  0 ) 

COMMON /P BOLD/ 3 10 .G 1 1 /DEL /HELP. DEL Z/ROLD/GSO . GS1 

C0MMON/RLK3/BS0.BS1/TIDY/EAST/ES/EKFK.EKEAST.TFK 

C0MMON/PLK3B/BI0.BI1 /TOY/N.NN.NNN/RLK2/X  I  .  WE /GGMS/GGMS 
COMMON /STORSCST/TSS IN ( 32. 21 >  .TSCOS ( 32.21 ) 

TPK  =  0.25  •  EKEAST 
PIOTWO  =  0.5  *  PI 
EPS1  =  0.0001 
DO  20  1=1*1 MAX 
R  =  PR  ( I ) 

IR  =  I 

00  20  J=  l.JMAX 
Z  =  PLZ(J) 

ARGP  =  h  -  Z  -  DEL 7 
IE(J  .EO.  1)  ARGP  =  0.0 
ARGM  =  H  -  Z 

5UMR=SUMI=SUMRR=SUMIB=SUMRV=SHMIV=SUMRBV=SUMIBV=0. 

DO  1  L  sl.N 

TRIGER  =  TSCOS (L . J )  *  RTO(IR.l.L) 

TRIGEI  =  TSSIN(L.J)  *  BTO(TR.I.L) 

SUMR  =  SUMR  -  TRIGER  *  BSl(l.L) 

SUMI  =  SUMI  ♦  TP  IGF  I  *  BSl(l.L) 

SUMRR  =  SUMRB  *  TRIGER  *  BIl(l.L) 

SUM  T  R  =  SUM  I B  -  TRIGEI  *  BIK1.L) 

SUMR V  =  SUMRV  -  TRIGER  *  BSO(I.L) 

SUM  I V  =  SUM  I V  -  TRIGEI  *  BSOd.L) 

SUMRR V  =  SUMRR V  -  TRIGER  *  RIO(l.L) 

SUM  1 RV  =  SUM  I R V  -  TRIGEI  *  BIO(l.L) 

1  CONTINUE 

SUMT(1)  =  SUMR  %  BUMT(l)  =  SUMRB 
SUM TV  < 1 )  =  SUMRV  %  RUMTV(l)  =  SUMRRV 
lE(U.EO.l)  GO  TO  54321 
FACT  =  ?./ (P»PI«ARGP) 

54321  CONTINUE 

FACT1  =  PIOTWO  /  AMINKROUT  -  R,  R  -  RIN) 

EEACT1  =  1.0  -  EXPEC-ARGP  *  EACH) 

ANSTR  =  RIN»SUMRR 
ANST  =  ROUT  *SUMR 
EL  =  0.0 
ICE  =  1 

2  ICE  =  ICE  ♦  1 

SUMM!T=SUMMTTR=VUMMIT=VUMM!TB=0. 
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no  375  M=1*N 

GMR  =  flTO  UR, ICE «M)  *  TSEXP <M, J, ICE) 

ARTIST  =  BSl(ICE.M)  *  GMR 
ARTISTR  =  RI1 (ICE,M)  *  GMR 
ARTISTV  =  GMR*RSO(ICE*M) 

VRTISTB  =  GMR*BI 0 ( ICE  *M> 

VUMMIT  =  VUMMIT  ♦  ARTISTV 
VUMMITR  =  VUMMITR  ♦  VRTISTR 
IE(J.EO.l)  GO  TO  93 
GO  TO  94 

93  GS  =  GTOUR.  ICF*M)*WF<M)  *TPK 
GSGS  =  GS*GS1 ( ICE«M) 

GSGSI  =  GS*G 1 1 (ICE.M) 

ARTIST  =  ARTIST  -  GSGS  %  APTISTB  =  ART ISTB  -  GSGSI 

94  SUMMIT  =  SUMMIT  ♦  ARTIST 
SUMMITR=  SUMMITB  -  ARTISTR 

375  CONTINUE 

mm  (ICE)  =  BUMT(  ICE-1)  ♦  SUMMITB 
SUMT(ICE)  =  SUMT (ICE-1 )  ♦  SUMMIT 
SUMTV ( ICE)  =  SUMTV ( ICE-1 )  ♦  VUMMIT 
BUM1 V ( ICE)  =  RUMTV ( ICE-1 )  ♦  VUMMITB 
EL  =  EL  ♦  EKEAST 
lE(J.EO.l)  GO  TO  95 
ELT  =  EL  ♦  E ACT  1 

EREC  =  AfiS (E ACT  *  EXPEI-ARGP  *  EL)  *  EEACT1  /  (EL  *  SUMT (ICE))) 
lEIEREC.LT. EPS)  GO  TO  R 
IEIICE.GE.JTOP)  GO  TO  fl 
GO  TO  ■?. 

95  ANSR  =  ROUT* SUMT (ICE)  ♦  .50 
ANSPR  =  R IN*BUMT ( ICE ) 

TESTER  =  ARS (1.0  -  ANST  /  ANSR) 

TESTERB  =  ARS (1.0  -  ANSTB  /  ANSRR) 

IEUCE.GE.JTOP)  GO  TO  100 

IEUFSTER.LT. FPS1.*N0.TFSTERR.L7.EPSI)  GO  TO  100 
ANST  =  ANSR 
ANSTR  =  ANSRB 
GO  TO  2 
100  CONTINUE 
NOEXI  =  0 

CALL  CEINT ( ARGM.  ROUT,  PR(IR),  EL,  CEM) 

StJMT  (ICF)  =  SUMT  (ICF)  -  CEM 
NOEXI  =  1 

CALL  CE?ERO(ARGM,  ROUT,  PR(IR),  EL,  CEM) 

NOEXI  =  0 

CALL  CE7ERO ( ARGP ,  ROUT,  PR(IR),  EL,  CEP) 

SUMTV (ICE)  =  SUMTV(ICF)  ♦  CEP  -  CEM 
CALL  CE INT ( ARGM,  RIN,  PR(IR),  EL,  CFM) 

BUMT(ICE)  =  BUMT(ICE)  ♦  CFM 
NOEXI  =  1 

CALL  CE7ERO ( ARGM,  RIN,  PR(IR),  EL,  CEM) 

NOEXI  =  0 

CALL  CE7ERO ( ARGP ,  RIN,  PR(IR),  EL,  CEP) 

BUMTVUCE)  =  RUMTV  (ICE)  ♦  CEP  -  CEM 
GO  TO  3 
B  NOEXI  =  0 

CALL  CEINT ( ARGP ,  ROUT,  PR(IR),  EL,  CEP) 
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SUMT  (ICF)  =  SUMT (ICF)  ♦  CFP 
NOEXI  =  1 

CALL  CF7ER0 ( ARGP .  ROUT*  PR(IR),  FL.  CFP) 
SIIMTV(ICE)  =  SUMTV(ICE)  ♦  CFP 
NOEXI  =  0 

CALL  CF INT ( ARGH «  ROUT,  PR(IR),  FL.  CFM) 
SUMT ( ICE )  =  SUMT(ICF)  -  CFM 
NOFXI  =  I 

CALL  CF7ERO ( ARGM ,  ROUT,  PR(IR),  FL.  CFM) 
SUMTV  ( ICE )  =  SUMTV (ICF)  -  CFM 
NOEXI  =  0 

CALL  CFTNT (AROP,  RIN,  PR(IR),  FL .  CFP) 
BUMT(ICF)  =  flUMT(ICF)  -  CFP 
NOEXI  =  1 

CALL  CF7ERO.ARGP,  RIN.  PR(IR),  FL.  CFP) 
BUMTV ( ICEi  =  BUMTVIICF)  ♦  CFP 
NOEXI  =  0 

CALL  CFINTIARGM,  RIN,  PR(IR),  FL.  CFM) 
BUMT(ICE)  =  BUMT(ICF)  ♦  CFM 
NOFXI  =  1 

CALL  CFZERO ( ARGM «  RIN,  PR(IR),  FL.  CFM) 
BUMTV ( ICF)  =  BUMTV (ICE)  -  CFM 
3  CONTINUE 

ANS  =  CMPLX (SUMT (ICF) »SUMI ) 

ANSB  =  CMPLX (BUMT( ICE) .SUMIB) 

VNS  =CMPLX (-SUMTV(ICE) .SUM  I  V ) 

VNSB  =CMPLX (-BUMTV (ICF) .SUMIBV) 

TSM(I,J)  =  ANS*R0UT  *  T IM ( I , J)  =  RIN* ANSB 
TSV(I.J)  =  VNS*ROUT  *  TIV(I.J)  =  RIN*VNSB 
IF(J.EQ.l)  TSM(I.J)  =  TSM ( I , J)  ♦  0.500 
?0  CONTINUE 
ENn 
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SUBROUTINE  CAL  TRIO 

C0MMON/RLK2/XI,WF/RLKl/H,R.FK,Pl/ST0RSCT!VTRSIN,TBC0S,TBEXP 
C0MM0N/T0Y/NQ01 .  ISYM , ICOR/T I DY/FAST 

DIMENSION  TBS  IN  (1?)  ,TRCOS(32)  .TRFXP(32,10>  .XI  (3?)  .XF<3?) 

EOUI VALFNCE  ( TRCOS ( 1 ) . TBEXP ( 1 )  ) 

COMMON/FS/FKFK.FKFAST.TFK 

COMMON/P  I T/JMAX  »  JMAXH* IMAX/OF.L/OFLR  »OELZ 

COMMON/STORSCST/STS  IN (32.21 > .STCOS (32.21 ) 

DIMENSION  STEMP ( 32)  ,CTEMP(32)  . ETEMP ( 32, 1 0 ) 

COMMON/S/RTO( 10.10.32) ,BT1 ( 1 1 ,  10.32)  ,QZ (32. ) 0) .SIGMA (32) .REST (128) 
COMMON/6/GTOU  0.1  0.32)  . STEXP  (  32. 20 . 1  0 ) 

EQUIVALENCE (STSIN(  1 )  ,TSSIN(1 ) )  ,  (STCOS ( 1 ) .TSCOS ( 1 ) ) 

EQUIVALENCE  (STEXP ( 1  ). TSEXP ( 1 )  ) 

OIMFNSION  TSSIN( 32.21) . TSCOS ( 32. 21 ) , TSEXP ( 32.20 •  1 0 ) 

FKFK  =  FK  *  FK  S  TFK  =  0.25  •  FK  J  FKF AST  =  2.0  *  FAST  *  FK 

TPK  =•  0.25  *  FKF  AST 

PTPKSQ  =  1.0  /  (TPK  *  TPO 

N  =  NQ01 

A  =  H  ♦  H 

00  300  L  =  1 .  N 

R  =  A  •  SIGMA (L) 

TBSTM(L)  =  SIN(B)  *  TFK 
TBCOS(L)  =  COS (R)  *  TFK 

300  CONTINUE 

A  =  -A.0  *  FAST  *  H 
B  =  FK  «  A 
TBMIJLT  =  EXP  (B) 

00  301  M  =  1 ,  N 
B  =  A  *  SIGMA (M ) 

TREXP (M.2)  =  EXP (B>  *  TPK 
00  302  ICE  =  3,  10 

TBEXP(M.ICF)  =  TREXP(M,ICF-I)  *  TBMULT 
302  CONTTNUF 

301  CONTINUF 
RETURN 
ENTRY  STTRIG 

HZ  =  0.5  *  OFLZ 
00  350  L  =  1  ,  N 
A  =  HZ  *  SIGMA (L ) 

STSIN (L  » 1>  =  S1N(A) 

STCOS  (L , i )  =  COS(A) 

PROO  =  STSIN(L.l)  *  STCOS (L  ,  1  ) 

STEMP (L )  =  PROO  ♦  PROD 

PROO  =  STCOS ( L . 1 )  *  STCOS (L ,  1  ) 

CTEMP(L)  =  PROO  ♦  PROO  -  1.0 
STSTN(L.l)  =  STSIN (L ,  1  )  *  TFK 
STCOS  (I.  «  1 )  =  STCOS  (L  .  1  )  •  TFK 
00  351  IR  =  2,  JMAX 

STSIN(l.IR)  =  STSIN(L.IR-l)  *  CTEMP(L)  ♦  STCOS  !L» IR-1 )  *  STEMP (L ) 
STCOS (L . IR)  =  STCOS (L « IR-1 )  *  CTEMP(L)  -  STSINIL ♦ IR-1 )  *  STEMP(L) 
3SI  CONTINUE 
350  CONTTNUF 
HZ  =  -HZ 

STMULT  =  EXP(HZ  «  FKFAST) 

00  352  M  =  )  ,  N 
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STFXP(M,1,1)  =  FXPCKZ  *  QZ  (M»  1  )  )•  *  TPK 
DO  353  ICE  =  2.  10 

STEXP(M.l.ICE)  =  STEXP(M.l. ICE-1)  *  STMULT 
ETEMP(M.ICE)  =  STEXP (M  » 1  *  ICE  >  *  STEXP <M. 1 ♦ ICE)  *  RTPKSO 
DO  354  IR  =  2*  JMAX 

STEXP(M,IR,ICE>  a  STEXP(M,TR-1,ICE)  *  ETEMP (M. ICE) 

354  CONTINUE 
353  CONTINUE 
352  CONTINUE 
RETURN 

ENTRY  SSTRIG 
DO  370  L  =  1  *  N 
DO  371  J  =  2»  JMAX 

STSINIL. J-l )  =  (STSIN(C.J-l) -STS IN (L i J) )  *  WF(L)  /  SIGMA(L) 
STCOS (L. J-l )  =  <STCOS(L. J-l )-STCOS(L»J) )  *  WE(L)  /  SIGMA (L) 

DO  372  ICE  =  2,  10 

STEXP(L»J-1«ICE)  =  (STEXP(L»J-1»ICE)-STEXP(L»J*ICE))  *  WF(L)  / 
*  QZ(L.ICE) 

372  CONTINUE 
371  CONTINUF 
370  CONTINUF 
RETURN 

ENTRY  TSTRIG 
JMAXP  =  JMAX  ♦  1 
DO  400  l  =  1»  N 
TSSIN (L.  1  >  =  0.0 

TSCOS  <L  *  1 )  a  TFK  *  WF(L)  /  SIGMA  (L) 

DO  401  J  =  2.  JMAXP 

TSSIN1L* J)  =  TSSIN (l. » J-l )  «  CTEMP(L)  ♦  TSCOS(LtJ-l)  *  STEMP(L) 
TSCOS(L.J)  =  TSCOS (l. * J-l )  *  CTEMP(L)  -  TSSIN(L*J-1)  *  STEMP(L) 
TSSINIL* J-l )  =  TSSIN(L.J-l)  -  TSS IN (L * J) 

TSCOS (L.J-1 )  =  TSCOS (L . J-l )  -  TSCOS(L.J) 

401  CONTINUE 
400  CONTINUE 

DO  402  M  =  1.  N 
DO  403  ICE  a  2,  io 

TSEXP (M. 1  *  ICE)  a  WF(M>  •  (1.0  -  ETEMP (M. ICE ) )  *  TPK  /  OZ(M.ICE) 
DO  404  J  =  2.  JMAX 

TSEXP (M»J* ICE)  a  TSEXP(M.J-l.ICE)  *  ETEMP (M* ICE) 

404  CONTINUE 
403  CONTINUE 

402  CONTINUE 
RETURN 
END 
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SUBROUTINE  CALBES 


DIMENSION  P7(20> .PL7(20) .PR(?0> .PLR(20) ,XI (3?) .WF (3?) 

OIMr^siON  BIO ( 10.3?) .BI1 (10. 3?>.GI0(10.3?).GI1 (10.32) 

DIM;" MS I ON  GS0(10.3?)«GS1 (10.32). RS0(10. 32). RS1(10«3?) 

COMk'ON/RLK  A/PR  .PLR.PZ.PLZ/P  IT/ JMAX  «  JMAXH. IMAX/PCC/ JTOP 
C0MMON/5/BT0 (10.10.3?)  .BT1  UI .10,3?) .07 (32. 1 0 > .S IGMA < 3? )  .REST < 12B ) 

COMMON /6/GTO (10.10,3?) .GT1 (11 .10,32) .RUSH (?8B0) 

COMMON/TOY /NODI »NOO?,NOD3/ROLD/GSO,GS1  /BLK1 /H» A.FK.P 1 /RLK2/X I »WF 
COMMON/RLK3B/BIO .P 1 1 /BROLD/GIO  »GI 1 /T IOY/F  AST/BLK3/RS0  »RS1 
COMMON /PAO/R IN. ROUT /GGMS/GGMS 
OIMFNSION  GGMS  < 1 0. 3?) 

N=N001 

FKFK  a  FK*FK  00015100 

FKFAST  =  2.*FAST*FK 
00  9  L=1 .N 

SIGMA(L)  =  X I (L )  ♦  0.5  •  FK 
OZ(L.l)  =  FAST  *  (XI (L )  ♦  X I (L )  -  FK) 

00  a  SIGMA(L)  *  SIGMA (L) 

0  =  SORT (FKFK  -  00)  00015700 

GGMS(l.L)  =  0 


T=0«R01JT  $  TI  =  D*RIN 
CALL  BESL(T.BSJO.RSJl.YO.YU) 

CALL  BESLUT.BIJO.BIJl  .YO.YU) 

BIO(l.L)  =  RIJO  *  Bll(l.L)  =  BIJ1  •  D 
BSO(l.L)  =  RSJO 
BS 1 (1 ,L>  =  BSJ1  *  0 
BTK1.1.L'  a  BIJ1  *  RIN  *  WF  (L >  /  0 
BTKIMAX  ♦  1.1. L)  =  BSJ1  *  ROUT  •  WF(L>  /  0 
9  CONTINUE 

0015  ICE«2* JTOP 
0016  Ha | «N 

OZ(M.ICF)  =  07 (M,  ICE  -  1)  ♦  FKF AST 
RHO  a  07 (M, ICE) 

OOR  a  RHO*RMO 

GMS  =  SORT (OOR  ♦  FKFK) 

GGMS(ICF.M)  =  GMS 
TM  a  GMS*ROMT  I  SM  =  PHO*ROUT 
TIM  a  GMS*R1N  $  SIM  a  RHO*PIN 
CALL  BESL (SM.GSTO.GSTI.Yl.Y?) 

CALL  BESL(TM,BST0,BST!  .YZ.YJ) 

CALL  BFSL(SIM,GITO.GIT1,YI,YO) 

CALL  BFSL(TIM.RITO.BITI.YK.YO) 

BIO(TCE.M)  a  BITO  %  RIl(ICE.M)  a  BIT1  *  GMS  $  GIO(ICE.M) 

Gil (ICE.M)  a  GIT  1 

BSO(TCE.M)  a  BSTO 

BS l (ICE.M)  a  BST1  *  GMS 

GSO(TCE.M)  a  GSTO 

GS1 ( ICE.M)  a  GST1 

BTK1.ICF.M)  a  PIN  *  BIT  1  *  WF  (M)  /  GMS 
RTKTMAX  ♦  1, ICE.M)  a  ROUT  *  BSTl  *  WF(M)  /  GMS 
GT1 ( 1 » ICE.M)  a  PIN  *  GIT1  *  WF (M )  /  RHO 
GTKTMAX  ♦  1. ICE. M)  a  POUT  *  GST1  *  WF(M)  /  RHO 
16  CONTINUE 
15  CONTINUE 

ODEI.R  =  (ROUT  -  PIN)  /IMAX 


00016000 


00016400 


00017100 

00017200 


00017700 

000171100 


GITO 


00018300 

00018500 

00018600 
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delr  =  .s*onFi.P 

RL  =  RTN  -  nnFLR 
no  S  IMP=1  <  TMAX 
RL  -  RL  ♦  OnELR 
no  i  l  ■)  *n 

D  a  OGMSCJj!  ) 

Tr  n*PR ( IMP ) 

T 1  a  n*(PR(JMP)  -  OELR) 

CALL  BESO(T.RTO<IMP,1,|> > 

IF(TMP.FO.l)  GO  TO  1 
CALL  RES1  (TNRTl  (IMP.l.L)  > 

RT1CTMP.1.L)  =  BT1  ( IMP*  1  »L  >  *  Rl  *  WF<L)  /  0 
1  CONTINUF 

00  S  ICF=2,JT0P 
00  ft  M=1 ,N 
RHO  =  07 (M  ,  ICE) 

GMS  =  GGMS < ICF*M) 

SM  a  PHO*PP ( IMP )  t  SMI  a  RHO* (PR ( IMP)  -  OFLR) 

TM  a  GMS*PR( IMP)  %  TM1  a  GMS*(PR(IMP)  -  DELR) 

CALL  RESO(TM,BTO(IMP»ICF.M) > 

CALL  RFSO'SMtGTO(IMP,ICE.M) ) 

IF ( IMP.FO. 1 )  GO  TO  ft 

CALL  BFS1(TM1,BT1(IMP,ICE»M)) 

BT1 (IMP. ICE. M)  a  BT1 ( IMP* ICE *M)  *  RL  *  WE(M)  /  GMS 
CALL  BEs;(SMl,GTI(IMP,ICE.M)) 

GT1 (IMP, ICE. M>  a  GTI (IMP, ICEtM)  *  Rl.  *  WF(M)  /  RHO 
6  CON T  TNI  IF 
5  CONTINUF 

DO  20  ICE  a  1,  JTOP 
DO  21  M  a  |,  m 
no  ??  IMP  =  1«  IMAX 

BT1 IIMP.ICE.M)  a  BTl ( IMP* 1  . 1CE.M)  -  RT1 (IMP* ICF.M) 
TF ( TCE  .EO.  1)  GO  TO  2? 

GTI IIMP.ICE.M)  a  GT U IMP* 1 ,  ICE.M )  -  GT 1 ( IMP, ICF.M) 
22  CONTTNIJF 
21  CONTINUE 
20  CONTTNIJF 
HFTURN 
FNO 


00016600 

00016700 


00018700 

00018800 


00018900 
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SUBROUTINE  CEINT (A*R«C*X* ANS) 

COMPLEX  ARG1.  ARG3.  El*  El 

COMMON/1CE/ICE.NOFXI/JOE/JTOP1/SAVE/ISAVF.FXMZR.CXMZI 
DIMENSION  A 1 RG (?)  *  A3RG<2).  EE1 (?)  *  EE3(?).  SAVER(?) 
EQUIVALENCE  (ARG1 * A1RG) *  (ARG3.A3R6).  (El. FED*  (E3.EE3) 
DATA  (PII  =  0. 15915494309) 

IE (ICE  .LT.  JT0P1 )  GO  TO  400 
EACT  =  ). 

IELAG  =  0 
NELAG  s  0 
GO  TO  1 
ENTRY  CEONE 

IE (ICE  .LT.  JTOP1 )  GO  TO  400 

EACT  =  1.  /  X 

IELAG  =  1 

NELAG  =  0 

GO  TO  1 

ENTRY  CEZERO 

IE ( ICE  .LT.  JTOP1)  GO  TO  400 
EACT  s  1.  /  X 
IELAG  =  1 
NELAG  =  1 

1  CONTTNUF 

EM  =  PII  /  SORT (R  *  C) 

A1RG(1)  =  A  *  X 

A1RG (2)  =  -  (B  -  C)  *  X 

A3RG ( 1 )  =  A1RG(1) 

A3RG (2)  =  -  (B  ♦  C)  *  X 
IE ( IELAG  .EO.  1)  GO  TO  2 
CALL  EXI(ARGl.El) 

NSAVF  *  ISAVE  S  SAVER(l)  =  EXMZP  «  SAVER(2)  =  EXMZI 
CALL  EXI ( ARG3.E3) 

GO  TO  3 

2  CALL  EX  12 ( ARG3*  E3) 

ISAVE  =  NSAVE  %  EXMZR  =  SAVER(l)  *  EXMZI  =  SAVER(2) 

CALL  EX  12 ( ARG1 «  El) 

3  CONTINUE 

IE(NFLAG  .EO.  0)  ANS  =  FM  *  (FF1 (2)  ♦  EF3( 1 ) ) 

IF(NELAG  .FQ.  1)  ANS  =  EM  *  (EE1  ( 1 )  ♦  EE3(2>> 

ANS  =  ANS  •  EACT 
RETURN 

400  ANS  =  0.0 
RETURN 
END 
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SUBROUTINE  EXK7.ANS) 

COMPLEX  Z«  ANS 
COMPLEX  ZZ.  AANS 
OIMFNSION  B (?) .  SUM (2) 

EQUIVALENCE  (ZZ.B).  (AANS.SUM) 
COMmoN/SAVE/ISAVE.EXMZR.EXMZI 
OIMFNSION  FACT (AO) 

DATA  (ISTORE  =  0).  (PIOTWO  =  1 .5707R63268) 

IE ( ISTOPE  .EQ.  1»  GO  TO  410 
DO  320  N  =  1.  40 

EACT(N)  =  -  FLOATIN  -  1>  /  ELOATIN  *  N) 

320  CONTINUE 
ISTOPE  =  1 
410  CONTINUE 
K  =  1 

GAMMA  =  0.S77215664R 

ZZ  =  z 

B1B1  =  B  ( 1 )  *  B  < 1 » 

B2B2  =  B (2)  *  R ( 2) 

II  =  B1R1  ♦  B2B2 
IEIU.GT. 100.0)  GO  TO  400 
ISAVE  =  0 

IE (B1B1  .GT.  B2B2)  SU  =  B<1>  ♦  0.5  *  B2B2  /  B(l) 

IE ( B 1 B 1  ,LE.  B2B2)  SU  =  ARS(  8(2)  ♦  0.5  *  B1B1  /  B(2)  • 
EN  =  6.0  ♦  3.5  *  SU 
NN  =  EN 

IE (NN  .GT.  40)  NN  =  40 

SUM ( | )  a  -R(I)  *  SUM (2)  =  -B(2) 

TERMR  =  -B ( 1 )  %  TEBMI  =  -R(2) 

00  300  N  =  ?«  NN 

X  a  (TEBMR*B(1>  -  TERMI*B«',>)  *  FACT  (N) 

Y  a  (TEPMR*B (?)  ♦  TERMI»R(l,»  *  FACT (N) 

TEPMP  =  X 

TEPMT  =  Y 

SUM ( 1 )  a  SUM(l)  ♦  TERMR 
SUM ( 2)  a  SUM (2)  ♦  TERM! 

300  CONTINUE 

ELNR7  a  0.5  *  ALOG(U) 

IE  (B  ( 1 )  .ME.  0.0)  FLNI7  =  AT  AN  (B  ( 2)  /  BID) 

IE (R  ( 1 )  .EQ.  0.0)  EINI7  =  SlGNCPIOTWOt  B(2)) 

SUM  ( 1  )  =  -GAMMA  -  EI.NR7  -  SUM  ( 1 ) 

SUM ( 2)  a  -ELNI7  -  SUM  (?) 

ANS  a  AANS 
RETURN 

400  CONTINUE 
ISAVE  =  1 

ONEOU  a  1.0  /  It 

SUM ( 1  )  a  1.0  S  SUM ( 2)  =  0.0 

TERMR  a  1.0  %  TERM!  a  0.0 

EEMMEE  a  0.0 

00  310  N  a  1*  4 

EEMMFE  =  EFMMEE  ♦  1.0 

FEE  a  EEMMEE  *  ONEOII 

X  a  -  EFE  *  (TERMR»B(1)  ♦  TERMI*R (2) ) 

Y  a  fee  *  ( TERMR*B ( 2 )  -  TERM !*B ( 1 )  ) 
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TERMP  =  X 
TERMI  =  Y 

SUM ( 1 )  =  SUM ( 1 )  *TE»MR 
SUM  <  ?)  =  SUM (2)  ♦  TERMI 
110  CONTINUE 

E  =  FXP(-B(1) 1 
X  =  E  *  COS (B (2) ) 

Y  =  -E  *  STN(R(?I) 

EXM7P  =  X 

FXM71  =  Y 

FACTO  =  (X  •  B(l)  ♦  Y  •  B(?))  *  ONEOU 
FACT!  =  (Y  «  B ( 1 )  -  X  *  R(?)>  *  ONEOU 
X  =  FACTR  *  SUM(l)  -  EACTI  *  SUM121 

Y  =  FACTR  *  SUM  (?)  *  EACTI  *  SUM  ( 1) 
SUM ( 1 )  =  X 

SUM ( ?)  =  Y 
ANS  =  A ANS 
RETURN 
END 
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SURROUTINE  EXI?(Z,ANS) 

COMMON/SAVE/ l SAVF»EXMZP»FXMZI/ ICE/ ICE,NOEX l 
COMPLEX  Z  »  ANS 
COMPLEX  ZZ,  AANS 
DIMENSION  R (?)  •  A<2) 

EO"»VALENCF  (Z7.B).  (AANS, A) 

AANS  =  ANS 

ZZ  =  z 
K  =  1 

TESTI  =  REAL (ZZ) 

TEST?  =  A IMAG  (7.7) 

IF ( TESTI  .EO.  0.0  .AND.  TEST?  .EO.  0.0)  GO  TO  1 
IE  (NOEXI  .EO.  0)  CALL  EXHZZ,  AANS) 

X  =  A(1)*B(1)  -  A (2) *B ( ? ) 

Y  =  A ( 1 ) *B ( 2 )  ♦  A ( 2 ) *R ( 1 ) 

IF ( I  SAVE  .FO.  1)  GO  TO  400 
E  =  EXP(-B(I ) ) 

EXM7R  =  E  *  COS (R ( 2 ) ) 

EXM7I  =  -E  *  S IN (R (?) ) 

400  CONTINUE 

A(l)  =  EXM7R  -  X 
A  (?)  =  FXM7T  -  Y 
ANS  =  AANS 
RETURN 
'  A ( 1 )  =  1.0 
A  (?)  =  0.0 
ANS  =  AANS 
RETURN 
END 
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5URR0UTINE  r,oc 

TYPE  ROUBLE  W.X  00000200 

COMMON /RLK 1/H»A»FK»PI /TOY /MODI »LL»MM/BLK2/XI ,  WF 
OIMFNSION  XI (32)  ,WF (32) ,W(31 ) ,X (31 ) 

OATA  <PJ  =  3.  141592653=9)  00000500 

OAT  A ( X  =  0. 1488743389816310000,  0.433395394) 29247D000, 

*  0.6794095682990240000,0.8650633666889850000, 

•0.97390652R5171 720000. 

*  0.076526521  1334973337550000. 0.2277858511416450780800000, 

*  0.3737060887154195606730000.0.5108670019508270980040000, 

*  0.6360536807265150254530000.0.7463319064601507926140000, 

*  0.8391 169718222188233950000,  0.912234428251 3259O586RDO0O. 

*  0.9639719272779137912680001,0.9931285991850949247860000, 

*  0 .99726386184948 156 j540000 »  0 .98561 1 51 1 54526 833540D000 . 


0.964762255587506430770000,  0.934906075937 7396891 70000,  00000700 

0.896321155766052123960000,  0.849367613732569970130000,  00000800 

0.7944837959679424  06960000,  0 . 7321 821 1 874028968038DOOO.  00000900 

0. 6630442669302 1520097D000,  0.567715757240762329040000,  00001000 

0.506899908932229390020000,  0.421351276130635345360000,  00001100 

0.331868602282127649770000,  0 .2392873622521 3707454D000,  00001200 

0.1 44 471961 58 279649 348D000,  0.48307665687738316234D-001 )  00001300 

OAT  A (W  =  0.2955242247147530000,  0.2692667193099960000, 

0. 2 1 90863625 15982D000,  0. 149451 349 15058 1D000, 

0.0666713443086880000, 


*  0. 1 52753387 1 307258506980000. 0.1 49 172986472603746788DOOO, 

*  0.1420961093183820513290000,0.131688638449176626  1980000, 

*  0.1181945319615184173120000.  0.1019301198172404350370000, 

*  0.0832767415767047487250000,0.0626720483341090635700000, 
*0.0406014298003869413310000,  0.0176140071391521 183120000, 

*  0.701861000947009660040-002.  0 . 1 6274394 T30905670605D-00 1 »  00001400 

*  0.253920653092620594550-001,  0.342738629130214331020-001.  00001500 

*  0.42835898022226680656D-001 ,  0.509980592623761762960-001,  00001600 

*0.586840934785355471450-001,  0.658222227763618468370-001 ,  00001700 

»  0.723457941088485062250-001,  0.781938957870703064710-001,  00001800 

*  0.833119242269467552220-001,  0.876520930044038! 1 1420-001 ,  00001900 

*  0.911738786957638847120-001.  0.938443990808045656390-001,  00002000 

*  0.956387200792748594190-001,  0.965400885147278005660-001)  00002100 

PFK  =  0.5  *  FK 

NNL  =  NODI/2  5  NRL  =  0 
IF (NOOl .EQ.20)  NRL  =  5 
IF (NODI ,EO . 32)  NRL  =  15 

00  111  1=1. NNL  00002400 

K  =  I  ♦  NBL 
X I ( I )  =  X(K)  *  PFK 
WF  ( I )  =  W  (K ) 

XI  ( I  ♦  NNL)  =  —  X  <  K »  •  PFK 
111  WF ( 1  ♦  NNL ) =  W ( K ) 

RFTIIPN 

FNTRY  TIME 

TIMFR  =  TIMELEFT  (K  ) 

PRINT  1.  TIMER 

1  FORMAT(*  TIME  LFFT  »F10.3*  SECONOS*//) 

RETURN 

FNO 
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SUBROUTINE  FARFLD 

DIMENSION  PTOP(IO) .FINS (20) .POUT (201 .ANS(bO) 

COMMON/LC/LCM AX .LCMAXH/U I ST/0 I ST/RAU/H IN* POUT /ANS/ ANS 
COMMON/HLKl/H.A.FK.PI/OEL/UELR.DELZ/BLKA/PR.PLR.PZ.PLZ 
COMMON/P  I T/ JMAX. JMAXH. IMAX/VELO/VEL/RRCC/RRCC.NPTS 
COMMON/ TOY /NOD  1 .ISYM.ICOR 

TYPE  COMPLEX  FF  I  .FF T .FF TB.RHOC. ADO  I T 1  .  SUM1 ,F F VT  .FF VS .BJ1 .BJOUT .CEP 

1  .CEM.P TOP. PI NS. POUT. ANS. 1  . SUM. AOD1 T . ARC I . VEL  .PHOT . V INS. VOUT . VTOP ♦ 

2  VHOT 

DIMENSION  VEL  (60) .PBOT ( 1 0 ) . V  I  NS ( 20 ) .VOUT (20 ) . VTOP ( 1 0 ) .SMUT  (10) 
DIMENSION  PZ  (20) .PLZ (20) *PR(20> .PLR(20) 

DAYA(ALICE=  B.6858B9b38) 

PRINT  700 
700  FORMAT ( 1H1 ) 

IF  (ICOR.EO.O)  GO  TO  41 
PRINT  40.DIST 

40  FORMAT  (//15X1 7HPATTtRN  AT  DIST  =  E 1 2 .4//20XSHANGLE 1 7X2HUB 
i  15X9HMAGNITUDE/) 

GO  TO  42 

41  PRINT  43 

43  FORMAT  (/13X28HFARFIEL0  PATTERN  AT  INF  INI T Y/20XSHANGLE 1 7X2HOB 

J  15X9HMAGNITUDE/) 

42  I  =  (0..1.)  $  RHOC  =  RRCC* I 

DD  =  l./(2.»NPTS)  $  UELTH  =  DD*PI  i  DELTHU  =  180. *DO 
NPT  =  NPTS-1  i  JMAXT  =  JMAXH 
IF (RIN.EO.O. )  JMAXH  =  0 

DO  2  J  =  1  *  IMAX 

Pt)OT(J)  =  ANSILCMAX  -  JMAXH  ♦  1  -  J) 

VBOT(J)  =  VEL (LCMAX  -  JMAXH  ♦  1  -  J) 

V TOP ( J )  =  VELIJMAXH  ♦  J) 

2  PTOP(J)  ■  ANS ( JMAXH* J ) 

DO  3  J=l. JMAXT 

POUT ( J )  =  ANS(IMAX.JMAXH.J) 

VOUT ( J ) =  VEL(IMAXUMAXHU) 

VOUT (JMAX  ♦ 1-J)  s  VEL (  LCMAX  *1  -JMAXH-IMAX  -  J) 

3  POUT (JMAX  *  1-J)  =  ANS (  LCMAX  *1  -JMAXH-IMAX  -  J) 

JMAXH  =  JMAXT 

IF (RIN.EO.O. )  GO  TO  966 
DO  965  J=). JMAXH 
PINS(J)  =  ANS (JMAXH*!— J) 

V1NS(J)  =  VEL ( JMAXH*1-J) 

V  INS ( JMAX*  1  -  J)  =  VEL (LCMAX  -JMAXH  ♦  J) 

965  P I NS ( JMAX ♦  1  -  J)  =  ANSILCMAX  -JMAXH  ♦  J) 

GO  TO  970 

966  00  969  J=1 .JMAX 
VINS(J)  =  (0.0. 0.0) 

969  PINS(J)  =  (0.0. 0.0) 

970  FF I =  ,5*FK*UELZ 

RO  =  FK*ROUT  $  HI  =  FK*RIN 

CALL  BESL(RO.BSO.BJO.YO.Yl)  S  CALL  BESL(RI.BSI.BJl.Yl.YO) 

BO  =  BSO*ROUT  S  B1  =  BSI*RIN 
H J I  =  RIN*BJ1  $  BJOUT  =  ROUT*BJO 
IF (ICOR.EO.O)  GO  TO  202 

HSI  =  RIN*RIN»HSI/DIST  t  BSO  =  ROUT*ROUT*BSO/DIST 
GO  TO  203 
202  BS I =BSO=0 . 
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203  BJOUT  =  BJOUT  ♦  l»HSU  4  BJI  =  BJI  ♦  I*BS1 
SUM  =  (0.0,0.01 
DO  37  K=  1  , JM AX 

AUDIT  =  <BJOUT»POUT (K)  -  BJI»PlNS(K>> 

ADO  1 T 1 =  KHOC* (VINS(K)*B1  ♦  VOUT(K)«BO) 

AOOI T  =  AUDI f  -  AUDI I 1 

37  SUM  =  SUM  ♦  ADO  I T 

FF1  =  Ft-  1  *SUM  4  FFVI  =  (0.0. 0.0) 

DO  38  K=1»IMAX 

RL  =  PH  ( K 1  -  .S«OELR  4  RU  =  HL  ♦  DELH 

TU  =  FK*HU  4  TL  =  FK«RL 

CALL  BES1 ( 1 U ,BU 1  4  CALL  BESl(TL.BL) 

38  FFVI  =  FFVT  -  (VTOP(K)  ♦  VBOT (K 1 1  * (HU*bU  -  RL*BL ) 

FFO  =  CABS(FFI  ♦  FFVT  *  «5*RH0C) 

902  DEG  =  0.0  4  FFOM  =  0.0 
PRINT  33.  DEG.  FFOM.  FFO 
DO  1  J=  1  »NPT5> 

THETA  =  J^DELTH 

IF (J.EQ.NPTS)  THETA  =  8999.*PI/18000. 

COSTH  =  COS ( THET A)  4  SINTH  =  SIN(THETA) 

ARGI  =  -I«FK*SINTH 
TT  =  RI*COSTH 

CALL  BESL(TT»BSl.hJl»YO,Yl) 

HI  =  R IN*BS 1 
TT  =  RO*COSTH 

CALL  BESL (TT.bS2.bJ2,Y0.Yl 1 

BJI  =  HIN*BJ1  4  BJOUT  =  R0UT«8J2  4  BO  =  R0UT*BS2 
IF ( ICOR.EO.O)  GO  TO  20S 

BSI  =  RIN*HIN«HS1/0IST  4  HSO  =  H0UT*R0UT*BS2/UIST 
BJOUT  =  BJ0UT*C0STH  ♦  I *0SO 
BJI  =  BJI«C0STH  ♦  I*BSI 
GO  TO  207 

205  BJOUT  =  BJOUT*COSTh 
BJI  =  B J 1 *COS  T  H 

207  FFI  =  .5*1*(CEXP(AhGI*0ELZ1  -  1.1/S1NTH 
SUM  =  (0.0,0.01 
DO  4  K=1,JMAX 
CEM  =  CEXP(ARGI*PL7(K)  1 

AUDIT  =  (BJOUT*POUT (K)  -  B J I*P INS (K 1 1 «CEM 
ADO  I T 1  =  ( VDUT (K 1 *BO  ♦  V  INS (K 1 *H I  1 *RHOC*CEM 

4  SUM  =  SUM  ♦  AOUIT  -AUUIT1 
FF I =  FF I *SUM 

SUM  =  (0.0,0.01 
FF  TB  =  ,5/COSTH 
ARGI  =  1*FK*H*SINTh 
CEM  =  CEXP(-ARGI) 

CEP  =  CONJG ( CEM 1 
DO  5  K  =  1 » I  MAX 

HU  =  PR  ( K 1  ♦  .5*DtLR  4  RL  =  PH(K)  -  .5*Dtl.H 
TT  =  FK«RU*COSTH  4  TTL  =  FK*RL»COSTH 
CALL  BESKTT.BU)  4  CALL  BES1  ( TTL.BL) 

RUBL  =  RU»dU  -  HL*BL 

AUDIT  =  -l»SINTH*(PBOT (K)*CEP  -  PIOP(K)*CEM) 

ADD  I T 1  =  HHOC*(VTOP(K)*CEM  ♦  VBOT(K)*CEP  1 

5  SUM  =  SUM  ♦  (AUDIT  -  AUDIT  1  I »RUHL 
FF  TH  =  SUM*FFTB 
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FFM  =  CAtlS<FFi  ♦  FFTH) 

DEG  =  J«UELTHU 
FFPHSUHt  -  FfF 
FFM  =  Al.  I  CE*LOGF  (FFM/FFO) 

PRINT  33.  DEG. FFM.  FFPRSURE 
33  FOPMAT <S> .F20 . 3.F20 . 3.E25. W > 

1  CONTINUE 

IF (OEG.GE .0. . ANO. ISYM.NE. 1)  GO  TO  900 
GO  TO  901 

900  UELTH  =  -UELTH  $  DEL  IHD  =  -OELThU 
GO  TO  902 

901  CONTINUE 
end 
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SUBROUTINE  S I M x I  MAT, MCT. RHS. ANSI 


DIMENSION  OHOLD 1 60 ) .MAT  (60*60) ,RHS (30) *  ANS (MCT  ) 

TYPE  COMPLEX  CC.CC2, OHOLD. RHS 

dimfns :on  c(?) *cx(P),cx2(?) 

TYPF  COMPLEX  MAT.ORIG, ANS,80,B2.R4.B6,B8,B10,B1 1 .B13.B15 
TYPE  INTEGER  V 

EQUIVALENCE  (R2.C) «  <CC*CX) « (CC?«CX?) 

DO  101  1=1. MCT 
00  101  J= 1 *MCT 
101  MAT ( I *MCT  * J)  =  MAT ( I « J) 

00  100  J=1.MCT 
OHOLO(J)  =  MAT ( J.MCT ♦ 1  ) 

100  MATU.MCT  ♦  1)  =  -  RHS(J) 

10  FORMAT (1X»C(E17.10,F17.10)  ) 

15  format (?sh  this  matrix  is  singular/) 

28  FORMAT ( 1 H 1 ) 

NCT=MCT ♦ 1 
JSING=JFIN=MCT 

IX  =  0S  R4=  (0.0*0. 0)1;  B1  1  =  (1  .0*0.0) 

JCT=MCT-1 
00  3  J=1  .  JCT 
KK=J*1 
GOTO  2S 

24  00  4  K=KK,MCT 

B8=MAT(K.J)/MAT(J,J) 

00  S  L=J.NCT 
R10=R8«MAT ( J.L) 

5  MAT(K*L)=MAT(K,L)-R10 
4  CONTINUF 

3  B1 1 =R1 1 *MAT ( J, J) 

B11=B11*MAT(MCT,MCT> 

LOWs-MCTt  M0=-1* 

00  6  INM=LOW*MO 
M=I ARS ( INM ) 

B0=-MAT (M.NCT) 

R2=MAT (M*M ) 

B4=(0. 0*0.0) 

IE (IX)  7.22.7 
22  IX=IX*1 
GOTO  8 

7  M02=-JFIN 

00  R  INN=LOw.M02 
N= I ARS ( INN ) 

9  B4=B4«MAT (M,N) *ANS (N) 

BO=RO-R4t  JF IN= JF IN-1 

8  IF (  C(l> *EO.O..ANO.C(2) .EO.O.)  13.29 

29  ANS(M)=R0/B2 

6  CONTINUF 
GO  TO  27 

25  V  =  J 

CC=MAT(J. J) 

IE (  CX ( 1 ) .EO.O.  .ANO.CX (2)  .EO.O.)  11.12 
11  IE(V.EO.JSING) 1 3, 14 
13  PR INT2R*  PRINT  1 5 
RETURN 


020 


040 

050 

091 

100 

120 

159 

160 
170 
180 
190 
200 
210 
220 
230 
240 
250 

00000255 

260 

261 

262 

270 

280 

290 

300 

310 

320 

330 

331 

332 
340 
350 

360 

361 
370 

410 

415 

420 

430 

440 

450 
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14 

V  =  V.l 

460 

CC2=MAT(V«J) 

465 

IE (0X2(1 ) .F  0.0. 

.ANO.CX2(2) .FO.O.l  11*16 

470 

16 

00  17  JJ=J.NCT 

480 

B6=MAT (J.JJ) 

490 

MAT (  I.JJ)=MAT(V 

.JJ) 

500 

17 

MAT (V* JJ) =86 

510 

Rll=-81 1 

520 

12 

JSING=JSING-1 

530 

GOTO  24 

540 

27 

00  200  J=1*MCT 

200 

MAT  ( J.MCT ♦  1 )  = 

OHOLD(J) 

00  102  1=1 «MCT 

00  102  J=1 «MCT 

102 

MAT ( T  *  J )  =  MAT< 

I.MCT.J) 

END 
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SURROUT  INF  RFS1  (X.H.)l  ) 

IF ( X.6T. 3. )G0  TO  1 
XT  =  <X*X>/9. 

fi.ll=  X*(.5  *XT* (-.56249985  ♦  XT* < .21093573  ♦  XT* (-.03954289  ♦ 

1  XT* ( .00*43319  ♦  XT* (-.00031761  ♦  XT*. 000011 09)>> )») 

RETURN 
1  XT  =  3./X 

FI  =  ,79788456  ♦ X T« (. 00000 1 56  ♦  XT* ( .01659667  ♦  XT* ( .000 1 71 05  ♦ 

1  XT* (-  .002495 1 1  ♦  XT*(. 001  13653  -  . 00020033*XT > ) )>) 

T 1  =  X  -  2.35619449  ♦  XT* (. 12499612  ♦  XT* (. 00005650  *XT* (-.00637879 
1  ♦  XT*(. 00074348  ♦  XT* ( .00079824  -. 000291 66*XT ))>)  ) 

SX  =  SORT ( X 1 

SX  =  l./SX 

RJl  =  SX*F 1 *COS ( T 1 ) 

RETURN 
ENTRY  RESO 
IE (X.GT.3.)  GO  TO  5 
XT  =  <X*X)/9. 

BJ1= ( 1 .  ♦  XT*<-2. 2499997  ♦  XT* ( 1 .2656208  ♦  XT* (  -.3163866  ♦ 

1  XT*(. 0444479  *  XT*(-  .0039444  ♦  XT*. 00021 00 ) > *>) ) 

RETURN 
5  XT  =  3./X 

FO  =  .79788456  *XT* (-.00000077  ♦  XT* (-.00552740  ♦  XT* (-. 0000951 2 
1  ♦  XT* ( . 00137237  ♦  XT* (-.00072805  ♦  XT*. 00014476) 111 ) 

TO  =  X-. 78579816  ♦  XT* (-.04166397  ♦  XT* (-.00003954  ♦  XT* ( .00262573 
1  ♦  XT* (-. 00054125  ♦  XT*(-. 00029333  <  .0001 3558*XT >>> > > 

SX  =  SORT ( X ) 

SX  =  l./SX 

RJ1  =  SX*E  0*COS ( TO ) 

RETURN 

END 
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SUBROUTINE  RESI.  (X«RJOtRJl*YO«Yl  ) 


1F(X.GT.3.)G0  TO  1 
XT  =  (X»X)/9. 


BJ0=(1.  ♦  XT* (-2.2499997  ♦  XT* ( 1 .2656208  ♦  XT*<  -.3163866  ♦ 

1  XT*(. 0444479  ♦  XT*<-  .0039444  ♦  XT*. 0002100) )))> ) 

BJ1=  X* ( .5  *XT* (-.56249985  ♦  XT* ( .21 093573  ♦  XT* (-.03954289  ♦ 
1  XT* ( . 00443319  ♦  XT* (-.00031761  ♦  XT*. 00001 1 09) ) ) )) > 


RETURN 
1  XT  =  3./X 

ro  =  .79788456  *XT* (-.00000077  ♦  XT* (-.00552740  ♦  XT 
1  ♦  XT* ( . 00137237  ♦  xY* (-.00072805  ♦  XT*. 00014476) ))> 
FI  =  .79788456  *XT* (. 00000156  ♦  XT* ( .01659667  ♦  XT* 

1  XT* (-  .00249511  ♦  XT*(. 00113653  -  ,00020033*XT) ))) 
TO  =  X-. 78539816  ♦  XT* (-.04166397  ♦  XT* (-. 00003954  ♦ 
1  ♦  XT*(-. 00054125  ♦  XT*(-. 00029333  ♦  .0001 3558*XT >) 

XT*(. 12499612  ♦  XT* (. 00005650  ♦ 
XT*(. 00079824  -.00029166*XT> ))) 


T 1  =  X  -  2.35619449  ♦ 
1  ♦  XT* ( .00074348  ♦ 


*(-.00009512 

) 

(.00017105  ♦ 

) 

XT*(. 00262573 
)  )  ) 

XT*(-. 00637879 
1 


sx  =  sort (x ) 

sx  =  l./SX 

BJO  =  SX*F  0*COS (TO  > 


BJ1  =  SX*F1 *COS ( T 1 ) 

RETURN 

ENO 
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FUNCTION  FN ( X ) 


OIMFNSION  AK<4)  .RK  (4)  .AE<4)  .BE  (4)  .A  (4)  ,P<4) 

OIMFNSION  FA (4) ,FB (4) ,GA<4) ,GBI4) 

OAT A  <FA=38. 027264.265. 187031. 135. 677320,38.102495) 

DATA  (FR =40 •021433*  32?. 62491 1 ,570.236280,157.105423) 

DATA  (GA  =  42.242855,302.757865,352.018498,21.821899) 

DATA  (GR=48.1969?7, 482. 485984,11 14.978885,449.6903261 
DATA (F2=-. 250000000 1  , IF4=0 . 0 1 041 6667 1  , (F6=-. 000231 48 1 1 
DATA(FR=0. 0000031001  , (F10=-. 0000000281  , <F3=-. 0555555561 
DATA (F5=0. 00 16666671, <F7=-. 0000283451  ,  (F9=0. 0000003061  , (Fll  = 

1-. 000000002) 

DATA  <AK=  0.09666344259,0.03590092383,0.03742563713,0.01451196212) 
DATA (BK=0. 12498593597,0.06880248576.0.03328355346,0.00441787012) 
DAT A (AE=0. 44325141463, 0.0626060 1220, 0.04757383546, 0.0 1736506451) 
DAT A (BE=0. 249983683 10, 0.09200 180037, 0.04069967526, 0.00526449639) 
FNTRY  Cl 
IFLAG  =  1 

1FIX.LT. 1.)  GO  TO  1 

27  XX=X*X 

FDEN  =  FB (4)  ♦  XX*(FB<3>  ♦  XX*IFB<2)  ♦  XX*(FB(1)  ♦  XX))) 

FNUM  =  FA  (4 1  ♦  XX* (FA (3)  ♦  XX*(FA(2)  ♦  XX*(FA(1)  ♦  XX))) 

F  =  FNIJM/  (X»F0EN) 

GNUM  =  G A (4 )  ♦  XX* (GA 1 3 )  ♦  XX*(GA(2)  ♦  XX»(GA(1)  ♦  XX))) 

GDEN  =  GB (4)  ♦  XX*(GB(3)  ♦  XX*(GB(2)  ♦  XX*(GB<1)  ♦  XX))) 

G  =  GNUM/(XX*GOEN) 

IF ( IFLAG. EO.O)  GO  TO  28 
FN  =  F*5 IN ( X )  -  G*C05 ( X ) 

RETURN 
1  XX  =  x*x 

CH  =  XX*  (F2  ♦  XX*  (F4  ♦  XX*(F6  ♦  XX*(F8  ♦  XXMF10  ))))) 

;AMMA  =  .577215664 
FN  a  GAMMA  ♦  LOGF(X)  ♦  CH 
RETURN 
FNTPY  SI 
IFLAG  =  0 

IFIX.LT.l.)  GO  TO  29 
GO  TO  27 

28  FN  =  -  F*COS (X)  -  G*SIN(X) 

RETURN 

29  XX  =  X*X 

PI  =  3.1415926536 

5H  =  X*(l.  ♦  XX*  (F3  ♦  XX*  (F5  ♦  XX*(F7  ♦  XX*(F9  ♦  XX*FU))>)) 

FN  =  SH  -  ,S*PI 

RETURN 

ENTRY  FLLIPK 

HOl.OA=l  .386294361  12  5  HOLOR  =  0.5 
DO  2  J=1 ,4 


A  (.1)  =  AMJ) 
2  R(.J)  =  BK  ( J) 
GO  TO  4 
ENTRY  ELLIPF 
HOLOA  =  1.0 
HOLOR  =  0.0 
00  3  J=1 ,4 
A ( J)  =  AE(J) 
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3  R(J»  a  RE  ( J) 

4  0  =  OM  =  1.  -  X 
FACT  =  I.OGFU./O) 

00  11  J=l*4 

HOLOA  =  HOLOA  ♦  A ( J) *0 
HOLOR  =  HOLOB  ♦  R(J)*G 
110=  0*0M 

FN  =  HOLDA  ♦  HOLOB*F ACT 
ENO 
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«IN=  0.040  H01JT  =  0.110  H=  0.100  FK  =  10.000  NU01=  32  IMAX  =  It  JMAA=  20 
THE  IMPEDANCE  COEFFICIENTS  E  OH  IH1S  PINO  TRANC  JUCEW  IN  UNITS  OF  WHO  C  A  AWE 


21 

= 

3.79389*000 

7.62228-001 

22 

= 

1.91 180-002 

2.36734-003 

23 

= 

0.14580-002 

2.62133-001 

2 1 P 

= 

3.04974-002 

-9.18393-002 

22P 

= 

3.36869-001 

-1.52465*000 

23P 

= 

5.36516-001 

-1.29543-001 

TIME  LEFT  151.173  SECONDS 
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WINb  THANsUUCEK 


»IN  = 


0.090 

POUT  =  0.110  M= 

NUM 

SP  HEAL 

1 

-1.68964067*007 

2 

-1 .67240098*007 

3 

-  1.63799449*007 

4 

-1 .5B65442b*007 

5 

-1.51816601*007 

6 

-1 .43285294*007 

7 

-1 .33023691*007 

8 

-1  .20902782*007 

9 

-1.06536050*007 

10 

-8.8494607 1*006 

1 

-5.94569027*006 

2 

-5.57339113*006 

13 

-5.21429150*006 

14 

-4.87234888*006 

15 

-4.54734618*006 

16 

-4.23676910*006 

1  7 

-3.93760985*006 

18 

-3.64749427*006 

19 

-3.3651 7701*006 

20 

-3.09074592*006 

21 

-1 .71773374*006 

22 

-1.09790283*006 

23 

-6.93093461 *005 

24 

-3.93304059*005 

25 

-I.62509b32*005 

26 

1.61794302*004 

27 

1.51970059*005 

28 

2.50266259*005 

29 

3. 14287642*005 

30 

3.45877678*005 

11 

3.45877678*005 

32 

1.14287642*005 

33 

2.50266259*005 

34 

1 .51970059*005 

35 

1.61794302*004 

36 

-1 .62509632*005 

37 

-3.93304059*005 

33 

-6.93093461*005 

39 

-1.09790283*006 

40 

-1. 71773374*006 

41 

-3.09074592*006 

42 

-3.36517701*006 

43 

-3.64749427*006 

44 

-3.93760985*006 

45 

-4.23676910*006 

46 

-4.54734618*006 

47 

-4.87234888*006 

48 

-5.21429150*006 

49 

-5.5733911 3*006 

50 

-5.94569027*006 

51 

-8.84946073*006 

52 

-1.06536050*007 

S3 

-1.20902/82*007 

54 

-1 .33023691*007 

55 

-1 .43285294*007 

56 

-1.51816601*007 

57 

-1 .58654426*00 7 

58 

-1.6  1/99449*007 

59 

-1 .67240098*007 

61) 

-  1 .68964067*007 

O.iOO  FK  =  10.000  N(JL»1=  32  JMAX  =  10  JMAX  =  20 


VEL  WEAL 


VEL  IMAblNAHY 


SP  1MAGINAWY 

-7.28632251*006 
-7. 1R043163*006 
-6.96905714*006 
-6.66287192*006 
-6. 2J 24 3279*006 
-5. 70743504*006 
-6.07518842*006 
-4.32697361*006 
-3.43745457*006 
-2.31425335*006 
-5.06481738*005 
-2.99199128*005 
-9.86962798*004 
9.32210914*004 
2.76379070*005 
4.51865036*005 
6.21313704*006 
7.86480965*005 
9.48761039*005 
1.10889949*006 
2.07624234*006 
2.53582577*006 
2.84441575*006 
3.07748136*006 
3.25955936*006 
3.40210606*006 
3.51133565*006 
3.59088113*006 
3.64290692*006 
3.66863508*006 
3.66863508*006 
3.64290592*006 
3.59088113*006 
3.51133565*006 
3.40210606*006 
3.25955936*006 
3.07748136*006 
2.84441575*006 
2.53582577*006 
2.07624234*006 
1  .  1  0889949*006 
9.48761039*005 
7.86480965*005 
6.21313704*005 
4.51865036*005 
2.76379070*005 
9.32210914*004 
-9.86962798*004 
-2.99199128*005 
-5.05481738*005 
-2.31425335*006 
-3.43745457*006 
-4.32697361*006 
-5.01518842*006 
-5.70743504*006 
-6.23243279*006 
-6.65287192*006 
-6.96905714*006 
-7. 18i)43163*006 
—7 .28632251 *006 


-1.03000000*000 
-1.03000000*000 
-1.03000000*000 
-1.03000000*000 
-1.03000000*000 
-1.03000000*000 
-1.03000000*000 
-1.03000000*000 
-1.03000000*000 
-1 .03000000*000 
-3.00000000-001 
-3.00000000-001 
-3.00000000-001 
-3.00000000-001 
-3.00000000-001 
-3.00000000-001 
-3.00000000-001 
-3.00000000-001 
-3.00000000-001 
-3.00000000-001 

9.70000000-001 

9.70000000-001 

9.70000000-001 

9.70000000-001 

9.70000000-001 

9.70000000-001 

9.70000000-001 

9.70000000-001 

9.70000000-001 

9.70000000-001 

9.70000000-001 

9. 70000000-001 

9.70000000-001 

9. 70000000-001 

9. 70000000-001 

9.70000000-001 

9.70000000-001 

9.70000000-001 

9. 70000000-001 

9.70000000-001 
-3.00000000-001 
-3.00000000-001 
-3.00000000-001 
-3.00000000-001 
-3.00000000-001 
-3.00000000-001 
-3.00000000-001 
-3.00000000-001 
-3.00000000-001 
-3.00000000-001 
-1 .03000000*000 
-1.03000000*000 
-1.03000000*000 
-1.03000000*000 
-1 .03000000*000 
-1.03000000*000 
-1 .03000000*000 
-1 .03000000*000 
-1 .03000000*000 
-1.03000000*000 


0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 

0.00000000*000 


THE  COMPLEX  WAUiAfloN  IMPEDANCE  IN  UNITS  OE  WHO  C  A  IS 
TIME  LEFT  147.382  SECONDS 


3.92364973*000  2.56547390*000) 
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FAkFIELU  HATTERN  AT  INFINITY 


ANGLE 

DB 

MAGNITUDE 

0.000 

0.000 

6.192189*005 

5.000 

-0.074 

6.139883*005 

10.000 

-0.296 

5.984582*005 

IS. 000 

-0.672 

5.731086*005 

20.000 

-1.210 

5.387252*005 

2S.OOO 

-1.921 

4.963764*005 

30.000 

-2.823 

4.473831*005 

3S.OOO 

-3.943 

3.932799*005 

40.000 

-5.316 

3.357707*005 

4S.000 

-6.997 

2.766838*005 

50.000 

-9.070 

2.179354*005 

5S.000 

-11.672 

1.615306*005 

60.000 

-15.033 

1.097012*005 

65.000 

-19.489 

6.‘-67744*004 

/  0.000 

-24.285 

3.780840*004 

75.000 

-23.670 

4.058214*004 

HO. 000 

-20.699 

5.713481*004 

85.000 

-18.968 

6.973370*004 

90.000 

-18.430 

7.419349*004 

TIME  LEFT  146.454  SECONDS 
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